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I. INTRODUCTION 


= ae ee eee eee eee 


The analysis of control systems, especially Bare 
systems, involves a great deal of salsulations. In the case 
of the classical analysis of single input / single output 
systems, these calculations. involve the manipulation of 
transfer functions. In the design process, wh2re the char- 
acteristics of the system componentS may change many times, 
the repetitious recalculation of transfer funstions can be 
vecy time consuning. Additionally, classical analysis tech- 
nigues are, generally, graphical in nature. The drawing and 
reirawing of plots for analysis can also take up a consider- 
able amount of the designer's time. The use of a computer 
to manipulate transfer functions and draw plots would allow 
the designer to spend a greater amount of his time 


designing. 


A. BACKGROUND 


The development of a computer aii for sortrol systen 
analysis is tertainly not anew idea. Many such aids 
already exist. However, most of these programs are special- 
ized in a particular type of analysis technique and many are 
not interactive. Thus, the designer must typicaliy itple- 
ment several of these programs during the design process. 
Each of these programs may have different methosis for input- 
ting data, different methods for controlling th3 program and 
different types o£ output. Several problems which exist 
with these programs are listed below. An attempt was made 
to minimize these weaknesses during the dsvelopment of 
DACSAP. 

1) The programs are not interactive. In these cases the 


user muSt, normally, type data ints a data file wsaeg 
a very specific format. 
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Zye Non-anteacactive Eee ues Viewed it EXeCcHEmon at data 
is net formatted preperly or will run to Completion 
BECGUGINg Medtargless OutpUe and will often leave the 
Operatoc with no idea where the error was made. 


3) Non-interactive aa do not offer the user any 
PegnalwcoOlmneouymullS, hepciine@e taltor the output to 
1is needs. ms a cCesult, eceams of tabular output may 
be produced when only a plot is desired 3c vice versa. 


4) The data file for oon-interactive programs must be 
re-accessed and edited and the program re-executec 
each time a change to the system iS made. 


5) Those programs which are interactive often require the 
user to kfow a specific language or set of colimands in 
order t)9 operate the program. 


6) In either 9f£ these types (of programs, the operator 
(Uist vem liens POSSeSSloOogy 2 User's manual whach 
outlines the input format. or the operating commands 
required to perform specific tasks. 


7) #=x‘The Bade a us are Specific to one type sf 2nalysis. 
The designer must, then, execute several Ereograms tc 
complete his. work, be ty OLUGeee neo mernput Gata and 
poral Papi eZLngG hinselr® swith the operating instruc- 
hrons LOE each. 

Although these programs may perform the reguired calcu- 
lations, a great deal of the user's time is Spent as a 


computer operator and less as a control system designer. 


B. PROBLEM OBJECTIVE 


The objective of this project was to design a comprehen- 
Sive program which would meet the needs of the designer 
without the problems listed in the previous section. To 
meet this goal, PiCenOolawe tg esequtLements were used as a 
guideline during the development of the program: 


1) The_ program should provide all of the most commonly 
used classical control system analysis techniques. 


2) The Shee Ree must be able to analyze both continuous 
MiGmcidnean sySsteices Thererere, transfer functions in 
the s, Z, w, wt domains must be handled. 


3) The program must be interactive. The user should be 
able to change any of the system parameters at any 
time without retyping all parameters. He should also 
be able to go from one analysis technigue to another 
Without having to re-enter any of the systen 
parameters. 


4) The Peg lae Should. be easy to operate or be "user 
yY a 


frien Errors input ky the user should be acy to 
detect and correct. The user should have complete 


ld 


control of the program, but not be required to know 
any special language or ccmmands. 


5) The ou a of the program should be selectable by the 
user and meet his immediate needs. 


6) The program should take advantage of state-of-the-art 
technology Sind ay eee computer power, memory and 
graphics capabilities. 


The reguirements listed akove are fairly general in 


nature. There are several ways of meeting these reqgquire- 
ments, particularly in the area of user friendly, interac- 
tive programming. The following chapter describes, in more 


detail, the methods used in meeting the objectives of this 
project. 


eZ 


II. DEVELOPMENT Or THE PROGKAS 
To meet the objective outlined in the previous chapter, 
the Digital and Analog Control System Analysis Program, or 
DACSAP, was developed. DACSAP is an interactive program 
written in FORTRAN to run on the IBM 370 model 3033 computer 
at the Naval Postgraduate Schocl. The FORTRAN programming 
language was chosen for the following reasons: 


feeen large library of FORTRAN Eunctions. and subroutines 
exist to perform routine numerical calculations. 


2) FORTRAN Bae can interface with the state-of-the- 
art graphicS package DISSFLA (see Appendix C). 


3) Since FORTRAN is a common programmin language, it 


will be easy for programmers to alter an improve 
DACSAP in the future. 


A. APPROACH TO THE PROBLEM 


A top-down programming technigue was used to develop 
DACSAP. First, the tasks that the program was to perform 
were defined. Then, each task was developed in greater 
detail until, finally, the very specific routines required 
to perform each task were implemented. Of course, a few new 
tasks were defined and routines created during the develop- 
ment of the program and were incorporated in parallel with 
existing routines. 

The tasks which the program performs may be broken into 
three basic categories: 


1) Those involved in manipulating the system's transfer 
PUNCtLOHS. 


2) Those involved in performing analysis of the system. 


3) Those involved in the control of the progran. 


The approach used in accomplishing each of tnese tasks is 


discussed in the following subsections. 
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1. System Manipulation 


DACSAP uses block @diageams to depict COnp@ael 
systems. Fach block in the diagram represents a component 
in the system and is described Ly a transfer function. 

Analysis is usually performed on the system's open 
or closed loop transfer functions. Therefore, DACSAP uses a 
single feedback control loop consisting of blocks aug 
forward and feedback paths for which the open and closed 
loop transfer functions may be calculated. Each Die@g@ awe 
transfer function is defined by a numerator and denominator 
polynomial in either coefficient or factored form. DACSAP 
makes avaiiable 20 blocks for transfer functions in the 
control loop. In most cases, this should allow the user to 
input the system exactly as it appears in the block diagram. 

It is common for a single input / single output 
system to consist of several nested feedback loops. These 
systems can be reduced to single loops by calculating the 
closed loop transfer function of the inner loops and 
including them aS separate blocks in the next outer loop. 
This method of block diagram manipulation is performed auto- 
matically by DACSAP and used to reduce multi-loop systems to 
a single loop for analysis. 

Routines were developed to allow the user to input 
multi-loop systems as they appear ina block diagram by 
Starting on the innermost loop and expanding to succeeding 
outer loops. DACSAP can be made to calculate the closed 
loop transfer function for the inner loop and place it 
anywhere in the next outer loop. 

It waS a primary objective to be able to change the 
system at any time with minimun effort. To this end, 
routines were developed to allow the user to change a single 
coefficient or root, redefine an entire polynomial, redefine 


an entire block or add new blocks. It was required that any 
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of these changes could be made without effecting any other 
part of the block diagram and that the resulting changes in 
mae system's open and closed loop transfer functions be 


automatically recalculated. 


2. Analysis Techniques 


One of the major objectives of the project was to 
make available, in a single program, as many control systen 
analysis techniques as possible. The following is a list of 
the anaiysis techniques chosen for inclusion in DACSAP; 

1) Root locus analysis 

2) Closed loop Bode frequency response analysis 
3) Open loop Bode freyguency response analysis 
4) Nyquist analysis 

5) Nichols analysis 

6) Time response analysis 

Each of these analysis techniyues may be perforned 
on digital or continuous systems. Thus, the user may choose 
to represent the system using Ss, 2Z, w or w' transfer func- 
tions. In all of the analysis technigues either the method 
of calculation or the way the data is presented differs 
depending on which domain is being used. So, routines were 
developed to handle analysis in any domain. This was done 
in all but one case. Frequency response of digital systems 
must be done in either the w or wt domains. 

It was also a requirement that the user be able to 
Switch from one analysis technigue to another without having 
to reenter any system parameters. This was easily accon- 
plished by using common block variables and arrays. During 
the development of the analysis routines, it was also found 
useful to be able to retain analySis parameters between 
executions of a routine. For example, if asystem is 
analyzed using Bode analysis, the user can leave the Bode 


routine, make changes to the system and rerun the Bode 
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analysis for the new system using the same freguency 
response parameters he had used before without having to 
reenter any of then. Of course, any or all of the parane- 


ters may be changed if the user so desires. 
3. Program Control 


Another of the project's objectives was to develop a 
user friendly method of program control which was flexible, 
but which did not require the user to know any special 
language or commands. This was accomplished using two 
concepts. First, the program is, to a certain extent, self 
Managing. This concert is applied only in those areas where 
total user controi is not felt to be necessary. FOr 
example, some programs require that the user instruct the 
program to calculate the open oor closed loop transfer func- 
tions after the system has been entered and prior to 
commencing any analysis. If the system is changed, instruc- 
tions must, again, be given to recalculate the transfer 
functions reguired for further analysis. 

DACSAP's operation is such that transfer functions 
are automatically calculated as they are needed to perforn 
analysis calculations. These transfer functions are also 
automaticaily updated when any changes are made to the 
system by the user. Although some control of the program is 
taken away from the user in these cases, the chance of 
making an error is also reduced and the user's time is spent 
more efficiently. 

The second concept used in program control is the 
use of option menus. A variety of menus are presented to 
the user throughout the program's execution showing the 
options available to him at that time. Self-management is 
also sometimes used in the ccntrol of the option menus. 
That is, the user is rarely given the opportunity to select 


which menu he will see next. Rather, the correct menu 1s 
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automatically presented to the user based on where execution 
is occurring in the program and which tasks have been 
completed to that point. 

The result of menu driven options is that the 
options available are presented in plain English so that the 
user does not have to know any special language or commands. 
Since the control of the menus is handied by the progran, 
the user is always presented with a menu of options when one 
is reguired and that menu only contains options which are of 
interest to him at that time. These concepts allow the 
program execution to flow with a minimum number of inputs 


mg tal fewer errors. 


B. ORGANIZATION OF THE PROGRAM 


As mentioned in the previous section, the top-down 
approach was used to design DACSAP. Subsequent to the 
design, Specific task reyguirements were defined. The 
routines to perform these specific tasks were written first. 
Thus, that the program was actually written from the bottom 
up. 

Low level routines were written based on their required 
Output. The input requirements of these routines, then, 
defined the output of the next higher routine. This process 
of building upward was continued until the required data was 
that data which must be supplied by the uSer. For this 
reason, the parameter input routines form the highest levels 
OF the program. 


The basic organization of the program is shown in Figure 


eee |< The lowest level routines are the analysis routines; 
these were written first. The calculation routines were 
written next and the parameter input routines last. The 


parameter input routines take inputs from the user and put 


them into forms which may be used by the caicalation 
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Figure 2.1 Program Organization. 

routines. The user controls the program primarily through 

these parameter input routines. Phe, these routines 


contain the instructions for presenting option menus and for 
changing and correcting parameters. 

The uvpermost controlling routine is the driving 
routine, called DACSAP, which provides the main option menu. 
This routine shifts program control to the parameter input 


routines based on the selection from the main menu. The 
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fasheda lines in Figure 2.1 show how control 1S Shifted. 
Once the required tasks have been performed by the analysis 
meta nes,  COoncrolwmws Shifted Eack to the driving progran. 
fifewe shitting of control to the caiculation and output 
routines is done autcmatically and is not controlled by the 
user. 

The flow of information is depicted by the solid lines 
mie ©rigure 2.1. The user supplies system data via the 
transfer function input routine and analysis parameters via 
the root locus, frequency response and time response paran- 
eter input routines. Information is taken from tnesa 
routines by the caiculation routinesS Which, in tura, provide 
miramrmation to the tabulation o£ plotting routines. 

The following chapter contains a more detaiied descrip- 


mon Of each of the branches Ssnown in TFigure 2.1. 
aa 


i 


III. PROGRAM DESCRIPTION 
DACSAP consists of a small driving program and 40 
subroutines. The driving program, which is called DACSAP, 
does little more than present the user with th2 main option 
menu (Fig 3.1) and call the appropriate subroutines based on 


the selection made from that menu. 





CS es ae ee ne 7 
| 
| | 
| DACSAP CPTIONS | | 
| | OPTION NO. | OPTION | 
1 INPUT TRANSFER FUNCTION (5) 1 | 

2 ROOT LOCUS ANALYSIS 

3 BODE ANALYSIS (OPEN L008 } | 

4 BODE ANALYSIS (CLOSED LOOP) 

5 | NYQUIST ANALYSIS i | 
| 6 NICHOLS ANALYSIS ba 
| 7 | TIME RESPONSE | | 

8 CHANGE BLOCK DIAGRAM 
l 9 | SAVE THIS PROBLEM el 
| 10 START A NEW PROBLEM lana 
1 HELP | a 
! 12 EXIT DACSAP | 

| 
| 
a a a 


Figure 3.1 DACSAP's Main Option Menu. 


All calculations, input/output functions and remaining 
program control functions are performed by DACSAP's subrou- 
tines. These subroutines fall into one of the following 
four basic categories: 

1) parameter input routines 
2) calculation and analysis routines 
3) plotting or tabulation (output) routines 


4) utility routines. 
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The four main branches of the program (transfer function 
input and three analysis Sequences) all contain subroutines 
from each of these categories, except for the transfer 
function input sequence which has no output routine. Each 
of the four branches will be described in greater detail 
later in this chapter. 

First, however, it 1S necessary to take a closer look at 
the control structure of the branches. Program self- 
Management is used to a greater extent when an analysis 


routine is first entered. Figure 3.2 depicts the relation- 
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Figure 3.2 Control Structure Before Parameter 
Initialization. 


ship between the main menu (DACSAP), the transfer function 


input routine and the analysis routines. Before analysis 
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can be performed, the system must be defined by using 
Option 1 of the main menu. Once the system has been input, 
the transfer function initialization flag is set. The 
significance of this is explained later. 

When an analysis technique is chosen from the main nenu, 
the user 1S presented with the option menu for that branch 
(coot locus, frequency response or time response). The user 
may select graphical or tabular output which will send hin 
to tne appropriate parameter input routine. Cnce all of 
these parameters have been entered, another Flag is set. 
Then, analysis calculations and output automatically soe@me 

After the first run through the analysis routines, the 
structure of control changes to that depicted in Figure 3.3. 
The change from the first structure (Fig 3.2) to this one is 
dependent on the condition of the transfer function and 
parameter initialization flags. Note that more direc 
control of the program is given to the user via the change 
menus. This allows the user to nake selected changes to any 
parameter and rerun the analysis in any sequence desired. 

The following sections briefly describe the function and 
capabilities of the four main branches. The basic mathemat- 
ical concepts used will be stated but not carried out in any 
detail. The reader is referred to the bdDibliography of this 
report for a list of texts which cover control system theory 


ind eta Li. 


A. TRANSFER FUNCTION MANIPULATIONS 


Transfer function manipulations are performed by several 
subroutines which permit the user to iiteractively input, 


Change, store or load a system's transfer functions. 
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Figure 3.3 Control Structure After Parameter 
Initialization. 


The purpose of the input routine, fTINPUT, is to 
permit the user to input a system's transfer functions just 
as they appear in a block diagran, then calculate the 
system's open loop transfer function for use by the analysis 
routines. Because of the large number of calculations 
required for transfer function manipulation, TINPUT is the 
largest of DACSAP's subroutines. 

The first time that the routine is entered, the user 
1S given the opportunity to input the transfer functions 
from the console or from a data file. Entry from the 


@onsole will be discussed first. 
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The system may consist of up to 20 components in the 
forward and feedback paths of a single feedback contioi 
loop. zach component must be defined by a transier function 
and by the path (forward or feedback) in which eames 
located. The numerator and denominator polynomials for each 
transfer function are entered separately and may be 
described by either polynomial coefficients or polynomiai 
roots. After each polynomial has been entered, the Yseruees 
presented with the coefficients or roots just entered and 
given the opportunity to change any which may have been 
entered incorrectly. 

Instead of typing in the transfer function pelwnee 
mials for a given block, the user may choose to load the 
block's transfer function from a data file (creating data 
files is discussed later in this section). If this opedena 
is chosen, the numerator and denominator polynomial coerfti- 
cients are pulled from the disk file, via the subroutine 
TLOAD, and placed in the appropriate arrays for tne current 
bLock. The purpose of this feature is to allow one toa 
design subsystems separately and store their transfer func- 
tions to disk. Then, when designing the outer loops of the 
systen, the subsystem may be loaded directly into its proper 
bihoc k. 

AS previously mentioned, when the transfer function 
input routine is first entered, the user may choose to load 
the entire system from a data file. Tf this optiongies 
chosen, all data pertaining to a stored control loop will be 
loaded from a user specified data file. This feature aliows 
the user to continue the analysis or design of a system that 
had been worked on during a previous terminal session 
without having to retype any system parameters. When the 
load is completed, the program will continue its operation 
as 1f the user had just finished entering the transfer func- 


t1OnS LEON themweonsole. 
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Once all transfer functions nave been entered, 
either from the console or a data file, andthe user is 
satisfied that they are correct, TINPUT calculates the 
steady state gain, roots and coefficients of the system's 
open loop transfer function. The results of these calcula- 
tions are stored in common variables and arrays to be used 


by the analysis routines. 


2. Qransfer Function Change Routine 


The instructions required to change or expand the 
meeek Gdlagcram are also contained in TINPUT. When any 
changes are made to the system, the routine autonaticaliy 
recalculates the parameters of the open loop transfer func- 
EA1N© Ill The option to change transfer functions may be 
selected from DACSAP's main menu or from any of the analysis 
routines' change menus. Those Chanyes which nay be 


pecformed on the block diagram are Shown in Figure 3.4. 





| 
| 
— te ee ce ee cee ce ee es Se eee Se Se es ee ee <a SED SD ses te ete me ete ee | 
1 CHANGE A taf. IN THe CURRENT LOOP | | 
ADD A NEW BLOCK TO THE CURRENT LOOP | | 
5 CHANGE DoF. iN AN INNER LOOP 
4 EXPAND TO AN OUTER LOOP { | 
| 2 NO CHANGES | 
6 Hee | 
— | 
Figure 3.4 Transfer Function Change Menu. 
If the first option is chosen, the user will be 


asked which block he wishes to change and is presented with 


the block change option shown in Figure 3.5. 
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Figure 3.5 Block Change Menu. 


Changing the block's numerator or denominator 
(options 1 or 3) means that the user simply wishes to alter 
one or more of the polynomial's coefficients or fcoots 
leaving ali other block parameters unchanged. 

By creating a new numerator or denominator (options 
2or 4), the user will be able to change the order of the 
_- polynomial and must enter all of the roots or coefficients 
for that polynomial. The form of the polynomial (factored 
or coefficient) must, however, remain the same as previously 
selected. 

By selecting the option to create a new numerator 
and denominator (option 5), all data concerning the current 
block is, essentially, erased and the user is sent to the 
transfer function input routine to re-enter a transfer func- 
tion for “that jalock= This permits one to enter transfer 
function polynomials of different form and order into any 
block of the control Loop. Just as one may initially enter 
transfer functions from a file, the user may choose to 
create the new numerator and denominator by loading then 


From a data filer 
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One may move a block from the forward path to the 


feedback path and vice versa (option 6). All other features 


about the block, including the block number, remain 
unchanged. Cniiyew tts focati1on watnin= =the Control loop is 
changed. 


TiPenecmseCCOUdeOpteLon fron Fagure’ 3.4 1S chosen, the 
user iS simply sent to the transfer function input routine 
to enter data for a new block. The block may be added to 
the forward or feedback path and may be entered from the 
console or from a data file. 

DACSAP was designed to operate on a single feedback 
eentrol Joop. However, since control loop transfer func- 
tions may be saved to and loaded from data fiies, multiple 
loop systems may be designed and analyzed in the following 
Manner. 

The user designs the inner loop(s) of the system 


Hicst, such as Loop 1 in Figure 3.6. Once the designer 15 





Figure 3.6 Sample Multiloop Systen. 


zy 


satisfied that Loop 1 meets his specifications, he saves the 
system parameters and starts a new problem by loading the 
inner loop's closed loop transfer function into an Wailer 
loop block and entering the remaining blocks as he wishes. 
Option 4 of Figure 3.4 simplifies this process. 

In option 4, those procedures are performed exactly 
as stated in the last paragraph except that the proper 
sequence of events occurs automatically and the user is 
prompted for each step throughout the procedure. The user 
may choose to place the inner loop transfer function into 
any block in the next loop. in the case shown in Figure 


3.6, Block 2 is selected and the outer loop appears as shown 





Figure: 3.7 Reduced Multiloop Systen. 


InN Fagure: ss. An important feature of this option is that 
the program keeps track of the block in the outer loop 
containing the inner loop transfer function. This becones 
important when this option (Option 4) is used in conjunction 
with Option 3. 
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In order to make changes in an inner loop, the user 
saves his outer loop, starts a new problem, loads the inner 
loop and makes the desired changes. Like Option 4, Option 3 
(Fig 3.4) automatically contrcls this sequence of events, 
prompting the user each step of the way. After this option 


has been chosen and all changes have been made to the inner 


loop, the user may return to the outer loop which now 
reflects the changes made in the inner loop. This is done 
by selecting Option 4 to re-expand to the outer loop. The 


user is asked if the outer loop should be entered from the 
console or a data file. Since the outer loop was saved 
under Option 3, its data file can now be loaded. In this 
case, however, since the program has kept track of the block 
into which the inner loop was previously placed, it automat- 
ically updates that block with the new closed loop transfer 
function from the inner loop. 

One can see that large multi-loop systems may be 
designed by creative use of Options 3 and 4 of the Transfer 
Function Change Menu. If each loop of the system is saved 
aS it is expanded, the user is free to return to any of the 
inner loops, make changes and have those changes automati- 
cally incorporated in the outer loops by successively using 


Option 4 to re-expand the systen. 


3. Transfer Function Save Routine 


The purpose of the Save routine is to store all 
transfer function data to a disk file under a user specified 
file name. This function is performed by the subroutine 
TSAVE. TSAVE is called either ty the user, through DACSAP‘s 
main menu, or by the transfer function change routine. 

The data file set up by TSAVE contains two sections; 
the block load section and the system load section. When 
TSAVE is executed, the user is asked if the closed loop 


meausrer function Shculd be calculated and saved. The 
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response to this question determines which transfer func- 
tion, open or closed loop, will be stored in the block load 
section of the data file. 

The user's decision to store one or the other is 
based on the type of system being saved. For example, if 
the designer uses open loop Bode analysis to design a 
cascade compensator, the open lcop (cascade) transfer func- 
tion is saved in the block load section. If the inner loeup 
of a multi-loop feedback control system had just been 
designed, the closed loop transfer function is saved so that 
it can be loaded into the outer loop. 

The system load section contains ali data wnich 
describes the entire control loop. These parameters include 
the domain of the transfer functions, the number of blocks, 
individual block transfer function polynomials, the open 
loop coefficients and roots and several other parameters. 
The system load section is always the same regardless of 
what is stored in the block load section. 

The resulting data file, then, contains ali ocMee 
data necessary to completely reconstruct the control ioop. 
The entire system may be recalled later for further analysis 
or just the open or closed locp transter function may be 


recalled for insertion into a block of another control loop. 


Be. ROOT LOCUS ANALYSIS 


The purpose of this routine is to calculate the complex 
roots of the system's characteristic polynomial over a user 
specified range of open loop gain values. The major subrou- 
tines involved in root locus analysis are DACRUT, which is 
the parameter input routine, and RTLCS, which calculates and 
plots pole, zero and root locations for the systen. The 
parameters which must be input by the user include a range 


of open loop gain values, the type of feedback (positive or 
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negative), the sample period for discrete systems and the 

dimensions and heading for the flot. 
PECincheNcntdrhymGentrOleesystem theory 1t 1S known that 

the characteristic equation for a single feedback control 


loop is 


1+ K N(s)/D(s) = 0 (ore) |) 
for negative feedback or 
1- K N(s)/D(s) = 0 (352) 
for positive feedback, where 
N(S) = numerator polynomial of the open loop transfer 
function 
D(s) = denominator polynomial of the open loop transfer 
function 


K = open loop gain 


Equations 3.1 and°3.2 may be rewritten, respectively, as 


i 
© 


D(s) + K N(s) (3.3) 


i 
-) 


Dis) - K N(s) (see) 


KILCS calculates an open loop gain increment based on 
the range of gain values and the plotting dimensions chosen 
by the user. This increment provides the greatest number of 
plotted points for each locus. For each gain value, kK, 
RTLCS calculates the coefficients of the characteristic 
Pelbynomial {Eqns 3.35 or 3.4), determines the roots of that 
polynomial by calling the IMSL routine ZPOLR, and plots the 


BOoOt locations. 
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Lines of constant damping ratio and natural frequency 
acre drawn on the root locus plct depending on which domain 
is being used. In the s-plane, lines of constant damding 
ratio are simply straight lines radiating from the origin 
and lines of constant natural frequency are circles about 
the origin. Because of this simplicity, s-plane plots are 
not cluttered by drawing these lines. 

On the w-plane, lines of constant damping ratio are 
curved lines radiating from the origin and lines of constant 
Natural frequency are ellipses about the origin. The 
w'-plane varies from one which looks like the w-plane to one 
Similar to the s-plane depending on the sampling Late = ieee 
enough lines are drawn on these plots to make tne user aware 
of the behavior of damping ratio and natural frequency for 
the system. The w'-plane roct locus ,plot in Figumeea 
Shows the lines: of constant damping ratio and natural 
frequency. 

On the z-plane, the stable region 1S within the unit 
circle and the lines of constant damping ratio and natural 
frequency follow a somewhat more complex pattern. The 
values associated with these lines vary with sample period. - 
These complexities make analyzing z-plane root locus piots 
more difficult. To make the task easier, the stable region 
and lines of constant natural frequency and damping ratio 
are drawn and labelled in detail for z-plane plots (Figure 
399) 2 

After the root locus plot has been drawn, the user may 


request detailed information about specific root locations. 


The program returns values for open 1oop gain, Natu 
frequency and damping ratio for these root locations. The 
user may, then, highlight the roots associated with a 


particular value of open loop gain. 


B2 
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C. FREQUENCY RESPONSE ANALYSIS 


These routines determine the frequency response of the 
system's open or closed-loop transfer function. Pier owe put 
of the frequency response routines may be dispiayed graphi- 
Sally or in tabular forn. Relative stability analysis may 
ke performed by examining the system's open loop frequency 
response using Bode flots, polar plots (Nyquist analysis) or 
dog magnitude-phase plots (Nichols chart analysis). The 
total system's steady state response to a Sinusoidal input 
may be analyzed using closed loop Bode plots. 

The major routines involved in frequency response anal- 
ysis are DACFRQ, the parameter input routine, and FREQR, the 
calculation routine. Three plotting routines, EPLOT, 
NYPLOT, and NICPLT, are used to graphically display the 
response. In the following subsections the input parameters 
which must be provided by the user are discussed and the 


calculations performed by the program are briefly described. 


1. Open Loop Analysis 


The systen's open loop frequency response is 
obtained by substituting jw for s in the open loop transfer 
function, GH({s) (in the discrete domains, wis replaced by 


jv or w' is replaced by jy'). The value of GH(s) at jw is 


GH(jw) = R(w) + j X(w) or) 


where 


R (w) 
X(w) = imaginary part of GH (jo) 


real part of GH (jw) 


GH(s) may also be expressed as 


GH(jW) = IGH(jw)i exp{jgd(jw)} (3.6) 
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where 


[GH (jo) 1 = RS a= are (3.7) 


and 


@(w) = tan-1{X(w)/R(w)} (3.8) 


are the magnitude and phase, respectively, of GH (ju). 

The coefficients of the open loop transfer function 
are known from MTINPUT. The user need only input, via 
DACFRQ, the open loop gain and the range of freguencies, w, 
over which the response is to be observed. 

For discrete systems, vor vy" are sometimes called 
"fictitious" frequencies and are not the same as real 


frequency, W. Instead, 


y= tan (wT/2) 
yt =§ (277 eee ange 17/2) 


If the user is working in a discrete domain, he may have the 
response represented as a function of either real or £icti- 
cious frequency. If real frequencies are desired, the user 
must supply the sampling period, T; the program wili then 
perform the frequency transformation. 

Once all of the above data has been input, FREQR 
evaluates GH{jw) at values of w within the frequency range 
selected by the user. Values c£ R(w), X(w), |{GH(w)f and 
@(w) are calculated and placed in common arrays. for use by 
the plotting or Eabu lation Eoutaines, 

Examples of Bode, polar and magnitude-phase plots 


are shown in Figures 3.10, 3.11 and 3.12, respectively. 
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Note that the Bode and polar flots are complete and that 
Bode and Nyquist analysis and design techniques may be 
applied directly to these plots, either in hardcopy or on 
the graphics monitor. The magnitude-phase plot, however, 
must be transferred to hardcopy and an appropriately sized 
Nichols chart (template) must be available to the user 
before Nichols chart analysis or design technigues are 
employed. 

For all of the technigues described above phase and 
gain margins and crossover frequencies are displayed when 
plotting 1s complete, if they can be determined over the 


specified frequency range. 


2. Closed Loop Analysis 


The calculations involved in closed loop frequency 
response analysis are very similar to those for open loop 
frequency response analysis, as discussed in the previous 
subsection. In this case, however, s is replaced by jwin 
the system's closed loop transfer function. The input 
routine, DACFRQ, prompts the user for the forward path gain, 
feedback path gain and the type of feedback. These parame- 
ters are then used by the subroutine CLOOP to calculate the 
coefficients of the closed locp transfer function polyno- 
Mials. The remainder of the parameters input via DACFRQ 
remain the same as described in the previous subsection. 
FREQR performs the frequency response calculations based on 
the closed loop transfer function and the user specified 
frequency range. Graphical presentation of the closed loop 


frequency response is in the form of a Bode plot. 


De. TIME RESPONSE ANALYSIS 


The routines described in this section provide the user 


with a time history of a closed loop system's response to a 
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user specified input. The user may choose to excite the 
system with either a unit impulse, a step function or a ramp 
mUAaCtLON. The method for calculating a continuous system's 
response is different than that for a sampled (discrete) 
systen. In either case, the user may display the response 
in tabular form or graphically as a plot of output Magnitude 


Versus time. 


1. Continuous Systems 


The major routines involved in evaluating continuous 
Peco time response are DACTIUM, the parameter input 
routine, STIME and DVERK, the calculation routines, and 
meEoT, the plotting routine. 

The parameters which are input by DACTIM are the 
type of input (excitation), the closed loop parameters and 
the amount of time over which the response is to be 
observed. The closed loop parameters are vaiues for the 
forward path, the feedback path gain and the type of feed- 
back (positive or negative). The subroutine CLOOP calcu- 
lates the coefficients of the closed loop transfer function 
based on these parameters and the open loop roots. 

In order to calculate the continuous time response, 
it is desirable to perform the inverse Laplace transform for 
the closed loop transfer functicn, thus obtaining an expres- 
Emem ©£or the outnut as a functicn of time. However, numer- 
ical methods for obtaining inverse Laplace transforms 
sometimes break down and cannot be relied upon in all cases. 
An alternate method is to generate a set of differential 
equations, and associated output equation, corresponding to 
the closed loop transfer functicn. A numerical metnod for 
solving that set of differential equations is then be 
employed and the output, asea s hunetron O& time, 1s 


calculated. 
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From elementary linear control theory, a transfer 
function, T(s), can be changed to an equivalent set of state 
eguations as follows. Given, for example, the transfer 


function 


Ts) =o eee ee (3.11) 
S) 7 26.0 Ee és 
s4 + a,s? + ans + a,S + a, 


- the corresponding matrix eguaticns are 











Pe) eee 
dp A) BO alt * | ol, Ga 
| Ke ~a,7 4a, ~da, ~a, | | X 4 | 1 | 
Yat) | Dow Dimitar! 
| X»5 | 
i (emt) 
se 
[| X, 
where 6 is the magnitude of the input at time, t. Eguation 


3.12 is the phase variable canonical form or companion 
Matrix form of the state equaticn and is based on the denon- 
inator coefficients cf the transfer function. Equation 3.13 
is the output equation and is kased on the zeroes' portion 
of the transfer function. 

The subroutine STIME calculates the coefficients of 
the matrices in Equations 3.12 and 3.13. These coefficients 
are then used by the IMSL routine DVERK to create and solve 
the corresponding set of differential equations for values 
of time over the user specified time interval. Common 


arrays are established for the output amplitude and time 
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values andare made available for the plotting routine, 
ie LOT . Figure 3.13 1s a sample continuous systen time 


response plot. 


Zz. Discrete Systems 
AS with continuous systems, it ais desirable to 
obtain an expression, as a function of time, for tke 
system's discrete transfer function. The inverse transform 


of a function written in the z-domain is obtained by exan- 
ining the function's Taylor series expansion. i eeecle magne — 


mone the z-transitorm, X{Z), Of a discrete signal, x*(t), is 


2 x (kT) z-k (3.14) 
(Oj hax tlyz=2 + x (27T)7—2 F25.¢ x (kIT) z-* +... 


X (Z) 


In the case at hand, 


x(Z) = T(z) 6 (2) (3.15) 


where T(z) 1s the system's closed Loop transfer function and 
6(z) is the z-transform of the input signal. TAZ LS 
obtained from the routine CLOOP. ZiUlbeealemlates 0(Z) for 
the type of input specified by the user and combines it with 
moze to form the rational function x (z). The function x (z) 
is expanded to an infinite power series by simple polynomial 
long division. For the convergent series, the coefficients 
mumtie  2-~ terms are the values of x (kT) in the pulse train 
mat). 
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In the case of systems described by w- or w'-domain 
eransirer functions, the same method is used. Pare st, 
however, the routine WINV is used to perform the algebraic 
inverse transform from w to zor w' to Z. Then the evalua- 
tion of the resulting z-domain transfer function is as 
described in the previous paragraph. 

The outputs of ZTIME are common arrays containing 
the values of x*(t) and t for the sampling intervals, kT. 
Since x*(t) is not a continuous function, values are plotted 
by TPLOT as discrete, unconnected points aS Shown in Figure 
3.14. 


E. HELP ROUTINE 


The word HELP appears as an option in the DACSAP main 
option menu, the transfer function change menu, and in each 
of the analysis routines' option menus. Li this Gptren is 
chosen from any of these menus, other than the DACSAP menu, 
several pages of information pertaining to the currently 
chosen routine appears on the screen. This information 
includes a brief description of the routine, a list of the 
required input data and very brief operating instructions 
for the routine. After the information has been read, the 
program returns the user to the menu from which HELP was 
selected. 

If HELP is selected from DACSAP's main menu, a brief 
description of DACSAP appears on the screen followed bya 
menu of options. From this menu the use may choose to read 
any of the previously mentioned HELP files or print a hard- 
copy of the entire HELP routine. A sample hardcopy printout 


is contained in Appendix B. 
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fen, OUMMARY 


In this chapter the capabilities of the major branches 
of DACSAP are covered as well as some of the concepts used 
in bringing about these capabilities. Appendix D contains a 
programmer's guide for DACSAP. An outline of each of 
DACSAP's subroutines is contained in this guide includins 
definitions of the main variables and arrays and notes on 
the specific functions of each subroutine. individuals who 
may wish to modify DACSAP in the future should find the 
programmer's guide and information in this chapter useful in 


analyzing the progran. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

The primary objective of this thesis was to develop a 
comprehensive computer aid for the design and analysis of 
linear control systems. The program was aimed at reducing 
or eliminating the large number of calculations involved in 
the design of control systems. To this end, the program was 
to allow for the entry of transfer functions as they appear 
in a block diagran, to incorporate all of the common Gonemen 
system analysis technigues and to permit the analysis of 


both continuous and digital control systems. 


Ae SUMMARY OF KESULTS 


The result of the thesis is the Digital and Analog 
Control System Analysis Program or DACSAP. DACSAP 1S a 
fully interactive design aid which incorporates a wide 
variety of methods for block diagram nanipulationgag 
control system analysis. The following list outlines the 


program's capabilities: 


1) Interactive, menu driven program control. 


Z) + BLOeK diagram manipulations including entry of 
transfer functions in coefficient or factored, ous 
open and closed One transfer function calculate 
and reduction of multi-locp systems. 


3) Continuous and discrete root locus, analysis including 
the anteractive examination of individual root loca- 
tions, tabular output and Gooe tocus pierce 


4) Continuous .and discrete eee response analysis 
including bilinear frequency transformations, open and 
closed loop Bode plots, polam picts” So anal- 
eae log magnitude-phase plots (Nichols analysis) 
and tabular outfut. 


5) Continuous and discrete time response analysis 
eta the response to Seep impulse and framp 
inputs, tabular OUtpueeand time Distory ever: 


6) Mass storage data. file routines which. allow the user 
to store he entire block dlagram pius the oOpenwan 
closed loop transfer functions for use during a 
Subsequent run of the program. 
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Pye GlLapnical OWUEpUt 1s produced using a state-of-the-art 

granaics naesage see appeni Ge wilechwcweates wtast, 

LGn ereSOruelon pilots a the user's terminal allowing 
real time analysis of the system and true interactive 
design. 

8) A HELP routine which provides the user with a, brief 
outline of operating instructions and available 
options for each of the program's major routines. 

The program is currently available for use by all users of 
the Naval Postgraduate School's IBM 370 mainframe computer 
and is Freing used by students and faculty for course work 


mad thesis research. 


Be. RECOMMENDATIONS 


During the development of DACSAP, an attempt was made to 
cover alarge number of the basic needs of the control 
system designer. There are Many more capabilities that 
could be added to the program which would enhance its 
utility in specific areas. Pnestollewging 1S a partial list 
of such capabilities: 

1) The ability to input sampling period as a variable in 
discrete system transfer functions, thus allowing the 
designer o observe the effect of varying sampling 
Cate without having to Becakewlate anys transfer 
functions. 

2) The pete ans OY to perform frequency response analysis 
of z-domain transfer functions. The routine. must 
include the automatic PEcionOhvdatlOle son | Z2Z-aonain 
transfer functions to the w-plane and the transforma- 
tion required to present the response as a function of 
real (unwarped) frequency. 


3) The se to sum transfer. functions and the ability 
to designate blocks in parallel feed forward paths. 


4) Routines ees an oriented toward the design of 
ain, lead and Bag comensators using root locus and 
requency respcnse techniques. 


5) Conversion of discrete transfer functions to contin- 
uous and vice versa. 


6) The ability to convert Syoucls represented in state- 
Space form to equivalent transfer functions. 


Unfortunately, the average user at the Naval 


Postgraduate School is allotted only one megabyte of virtual 
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memory on the IBM 370 and DACSAP (in combination with the 
DISSPLA routines) requires nearly all of that memory to run. 
Therefore, any additions to the program would have to be at 
the expense of some existing routines or the structure of 
the program would have to be changed. 

Since the list above involves mostly transfer function 
Manipulations, one possibility for restructuring the program 
wouid be to separate the transfer function manipulation 
routines from the rest of the program (analysis and plotting 
routines). The results of the transfer function routines 
could te placed in a data file which could then be read by 
the analysis routines. However, this would reduce the flex- 
ibility and ease of operation of the program and increase 
its execution time. 

Thus, room does exist for increasing the capabilities of 
DACSAP. However, some trade-offs will probably have to be 
made if additions are to be made to the program, and that 
program is to be availabe to all users at the Naval 


Postgraduate School. 


ai0 


eeiiammcitemetiitimntinametittendtiminemest 


USER'S MANUAL 


Digital and Analog Control System Analysis Progran 
DACS AP 


DACSAP is an interactive aid for the design and analysis 
©: linear, slugresinput 7 single “output, (SISO) COne ro. 
systems. The program is designed to manipulate transfer 
functions from a control system block diajram and create a 
variety of output plots or tabulated data suitable for anal- 
ysis. Analysis technigues which may be used include root 
locus, Bode, Nyquist, Nichols and time response. 

The program is fully interactive and easy to use, 
requiring the user to know no special language or ccmmands. 
All data entry is precmpted and options are presented to the 
user via a set of option menus. For these reasons, this 
Manual is meant to be more of a tutorial for first time 
users. It is hoped that the user finds the program so easy 
to use that he no longer needs to refer to this manual after 
the program has been run a time or two. Towards this end, 
HELP routines may be accessed during program execution. The 
HELP routines briefly summarize the major points contained 
in this manual and may be used to refresh the experienced 


user's memory as to the operating procedures for DACSAP. 
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Aj. GETTING STARTED 


At the U.S. Naval Postgraduate School, DACSAP resides as 
a DISSPLA text module ona separate disk which may be 
accessed by typing the word CONTROLS after logging on. 
Since DACSAP, in combination with the DISSPLA subroutines, 
takes nearly one megabyte of memory to operate, the user 


should define his virtual memory as one megabyte (1024 K) 


prior to accessing the CONTROLS disk. If this is not done, 
the CONTROLS disk MANAGE EXEC will define the proper 
storage, then halt execution. The user will then be 


reguired to reenter the CMS oferating system and reaccess 
the CCNTROLS disk. Once on the CONTROLS disk, DACSAP is 
executed by selecting it from the MANAGE EXEC menu. 

When program execution has Ltegun, the user will be asked 


to select an output device as shown below. 


GRAPHICAL OUTPUT DEVICE OPTIONS: 


lebih Olle 
2- TEK 4113 OR 4114 
Se, LBteesz 73 


fe DESSPeA NEPARIECE 
PLEASE SELECT A GRAPHICAL OUTPUT DEVICs (1-4). 


Since DACSAP is meant to be an interactive design tool, it 
is best utiiized when output is to a high speed, high reso- 
lution device such as the TEK 618 or TEK 4113/4114 witha 
high speed moden. If one wishes to input a large systen, 
make changes to a system or oktain tabular data prior to 
performing graphical analysis, and do so without tying up a 
graphics terminal, then the IBM 3278 terminal mav be 
selected. Graphical output to these terminals is possible, 
but is in the form of an unlabelled, low resolution, line 
printer style plot. The general shape of curves can be 
observed, but detailed analysis of any kind is practically 


impossible. 


2 


Output to a DISSPLA metafile is available to those users 
who desire hnardcopy, RUD bEcattgoneGudlt ty Solotseiomasuse in 
theses, technical reports, etc. After the metafile has been 
formed tne user is free to select any of a number of output 
devices available through the DISSPOP routine. The actual 
plots produced are of the proper size for thesis figures and 
would, therefore, fit comfortably in any document written on 
ao) x iit inch paper. The figures for this document were 
produced using this option. The program may be run from any 
terminal if this option is chosen since no graphical output 


comes to the terminal. 
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B. TRANSFER FUNCTION INPOT 


Before any analysis may be performed, it is first neces- 
Sacy to define the systen. Each loop in the system nay 
contain up to 20 blocks each of which may be of 20th order 
or less. However, the entire system must be of 100th order 
or less. In order to define the systen, select Option 1 
from the DACSAP Option menu belcw: 


ee ee ee ee ee ee eee ee ee eee ee ee ee ee ee ee ee ees ee ee ce ee ee ee ee ee 


| “OPil Oley o. | OL PLON | 
> | 1 i INPUT TRANSFER FUNCTION (TS) i < 
| 2 i ROOT LOCUS ANALYSIS i 
5 BODE ANALYSIS (OPEN LOOP) 
vl BODE ANALYSIS (€ LOSEDREOOr) 
5 NYQUIST ANALYSIS 
6 NICHOLS ANALYSIS i 
ii TENE RESPONSES 
8 CHANGE BLOCK DIAGRAM 
Q SAVE THs PROBE { 
10 START A NEW PROBLEM | 
11 HELP 
| 12 | EXIT DACSAP i 


Selecting Option 1 will preduce the following menu of 


transfer function input optione. 


<> OD oD ee ee ee ee ee ee ee ee eee ee ee eee ee ee eee ee ee eee eee ee a eee Se ee oe ee 


| OPTION NO. | OPTION | 
i 1 | ENTER XFER ES Cay FROM CONSOL2 | 
2 { LOAD TRANSFER FUNCTION(S) FROM A FILE | 
3 RETURN TO DACSAP MENU i 
4G Pe i 


Option 2 from the menu should be selected if the user 
has saved a system during a previous terminal session and 
would now like to continue work on that system. The user 
will be prompted for the filenane, then the system will be 
loaded and transfer function infut will be complete. 

Option 1 should be selected if the system is being 
defined for the first time and will result in the following 


questions: 
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Pieroni leH DONATI SARE YOU WORKING? {35,2,9%) OR WP) 
WHEeH LOOP ARE YCUSPREPARING TO ENPUT? 
MOw MANY BLOCKS ARE IN YOUR OPEN-LOOP? 


It should be noted that blocks which consist of purely a 
gain should not be entered as transfer functions. Tne 
program requests values for open loop gain, forward path 
gain and feedback path gain whenever they are needed for 
calculations. The following questions are asked for each of 


the blocks in the open loop: 
TRANSFER FUNCTION NO. 1 


moeo Hts BLOCK IN THE FORWARD OR PEEDBACKMPAT I? (F OR 3B) 


DO YOU WISH TO INPUT THE TRANSFER FUNCTION FOR THIS 
BLOCK FROM A DATA FILE OR FROM THE CONSOLE? (D OR C) 


If entering the block's transfer function from the console, 


the following questions will be asked: 


WHAT IS THE ORDER OF THE NUMERATOR POLYNOMIAL? 
fone Loe thE ORDER OF THE DENOMINATOR? 


(E eee) ee PONG ETON er Reon so SOR CORSET LeCLENT FORM? 


ie Leanster Functions in BRactored Form 


The numerator and denominator of each block must be 
of the same form. For transfer functions in factored fora, 


input is prompted as follows: 


WHAT IS THE NUMERATOR GAIN CONSTANT? 
WHAT ARE THE ROOTS OF THE NUMERATOR? 


ROOT NO. 1 
REAL PART = 


IMAGINARY PART = 


The last three lines are repeated a number of times equal to 


the order of the polynomial. After the numerator roots are 
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entered, they are displayed as follows allowing the user to 


check for and correct any typing errors. 


> ee ee es eee es eee ee es es ee ee eee ee eee es ee ee ees ee ee eee ee 


> <a ce ee ee ee ee ee eee eee eee eee eee eee eee eee eee ee 


| LOOP NUMBER 1 | 
BLOCK NUMBER 1 | 


| ITEM NO. | NUMERATOR ROOTS 
| 1 -0. 0300000 0.0000000 3 | 
2 -4. 3000000 2.5780000 Gam 
3 -4. 3000000 -3.5780000 J. | 
| Uy CONSTANT = 25.0000000 | 


ANY CHANGES TO THESE PARAMETERS? 


Entry of the denominator polynomial is identical to 


that of the numerator. 


2. Transfer Functions in Coefficient Form 


If the transfer function is in coefficient forn, the 


user will be prompted as follows: 


NUMERATOR COEFFICIENTS. 
COEFFICIENT SCr Ss 
COETP REI SNT .CF ssa 
COEFF IECSENT Cr SS) aes 


Oo —- NA) 
ll 


Naturally, the number of coefficients reguested will corre- 
spond to the order of the polynomial (the case above would 
be of second order). After they have been input, the coef- 
ficients are displayed as follows and may be corrected if 


necessary. 
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—_ an Pan aap ame ewanm ae SEP ewamanm am am aR ee ewe aman SE aw SRF aan SP aww 2 RP aD Ee ewan Bap ewe ame ap 4 SRP aD ewe ew em aw ew es aR am a a 


| LOOP NUMEER 1 | 
BLOCK VNUMB ER a2 | 


{| ITEM NO. | NUN RALOR = COn EP aeLeN TS | 
| 1 le UOOUO 00M S ** 2 | 
| Z SOeseououg, 5 ** | 

3 | He oCuoCg > Ss =x 0 
mo PCHANGE S STORTHESE PAR ARETIERS ? 


Entry of the denominator polynomial is identical to 


that of the numerator. 


3. Transfer Functions ina Data File 


= Se ee —_— oe SS ca os eee ee = > ae oo ae oe 


If an inner feedback loop, a compensator, a filter 
or any other subsystem had been designed and saved using 
DACSAP in a previous terminal session, its transfer function 
could now be loaded into the appropriate block of the 
current system. This allows the user to design subsystems 
independently of the overall system then install them later. 
The user iS Simply prompted for the filename under whica the 


subsystem is saved. 


ay 


C. ROOT LOCUS ANALYSIS 


Root locus plots ora tabulation of root locations nav 
be obtained by selecting Option 2 from the DACSAP option 


menu shown below. 


> <a ce ee ee eee ee ee eee eee ee See eee eee eee eee eee eee eee eee eee eee 
ED SG DO Se ee ee ee ee ee ee ee ce ee eee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee 


FER _PUNCTION (S) 


{ODEN LOOP) 
CLOSED LOOP) 





NS 
S ANA 
YSIS 
ALYSIS 
NAYS 
NALYS re 
ONSE 
OCK DIAGRAM 


This selection will yield the following menu of root locus 


options: 
| ROOT-LOCUS OPTIONS j 
|} OPTION NO. |] OPTION | 
| 1 | ROOD LOCUS@: LO. | 
| Z | TABULATED DATA | 
3 RETURN TO DACSAP MENU 
4 EXIT TLCACSAP 
5 | HELP 


If a plot is desired, Option 1 should be selected. The 


user is then prompted for information in three areas as 
follows: 


GAIN PARAMETERS 
MIN GAIN VALUE 
MAX GAIN VALUE 
POSITIVE OR NEGATIVE FEEDBACK? (P OR N) 
TIC MARKS? (Y OR N) 


WHAT IS YOUR SAMPLE PERIOD? 
a (discrete systems only) 


SO 


PLOTTING DIMENSIONS 
X MAX 
X MIN = 
Y MAX = 
Y MIN = 


PLOT HEADING 
HOW MANY LINES OF HEADING WOULD YOU LIKE? (4 MAX) 


A MAXIMUM OF 32 CHARACTERS PER LINE IS ALLOWED. 
LINE 1 = (cepeated as required) 


If the user is unsure of how to dimension the root locus 
plot, it may help to first select Option 2 of the Root Locus 
Options menu. This will produce the following information 


at the terminal; (Sample values) 


] DEGREE | NUMERATOR | D=ZNOMINATOR | 

i ub | | 1.0000000 | 
3 -9400000 | 
y Zoe CoO UUe Ceo ©0002 
1 25.7499 847 One 07 9997 

i 0 i Oe 7499539 | Oe 575999 7-7 


POLE-ZERO LOCATIONS 
j2 ON Raps 


Ro ADReEART IMAG. PART 

eco 99 24an+ 01 0.00000000E+00 

-~0.40000027E+00 0.00000000E+00 

OR ooo Det 00 ORES 103 1) A See ONO, 

Ore 2 SISIS5E + 00 Ms 572012 9E+00 

ALL ZEROES ARE AT ECNPUNGET WecxCeP merOR THE FOLLOWING 

REAL PART IMAG. PAKT 

=O 799999530 E+00 0. 00000000E+00 

-0.30000009E-01 0. 00000000E+00 


Toupee LO) Lip AstABULAR SOUIPUT OF ROOT LOCATIONS? Gr SOR en ) 


A detailed iisting of root locations may be obtained by 
answering yes (Y) to the last question. By answering no (N) 
to the question, the user is returned to the Root Locus 


Options menu and may choose the option to produce a plot 


ao 


(Opt reonst)=. The takulation of pole-zero locations Shemma 


help in choosing the dimensions for the plot. It should be 
noted that the plot area is always square. So, the dimen- 
sions should also be chosen © pe syuare 


(i-e., XMAX-XMIN=YMAX-YMIN) af true angular representation 
is desired. 
When the root locus plot has been completed, the 


following question will appear on the screen: 


WOULD YOU LIKE INFORMATION ABOUT A PARTICULAR POINT 
ON THE ROOT-LOCUS? (Y OR N) 


By replying yes (¥), informaticn concerning any root loca- 
tion may be presented by responding to the following 
PEOIP tS: 
WHICH ROOT LOCATICN? 
REAL PART = 
IMAGINARY PART = 


The following information will be provided for that root 


location: (sample values) 


GAIN = 1.5 746107 
NATURAL FREQUENCY = a5) Odes RADZS2 
DAMPING RATIO = 0.7071068 


ANOTHER POINT? (YOR N) 


More root locations may be examined by responding with a yes 
(Y) to the last question. When a negative cesponse is 


input, the following question is asked: 


WOULD YOU LIKE TO HIGHLIGHT A PARTICULAR SAIN VALUE? 


If a yes (Y) is input, the user is asked for a gain value 
and a small box is drawn around the root locations associ- 
ated with that gain value. At this point the root locus 
plot is complete. Pcle-zero locations and the stable region 


crossover gain are then presented at the terminai. The 
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Stable region crossover gain is that jain at which the root 
locus first crosses from the stable region to the unstakle 
region. 

Lines of constant damping ratio and natural frequency 
are drawn on plots for discrete systems. BOG ijw- and 
w'-domain plots only a few lines are drawn, just to make the 
user aware of the behavior of these characteristics. Due to 
the complexity of these lines on the z-plane, they are drawn 
and labelled in greater detail. This should help make these 
plots easier analyze. Of course, the most detailed informa- 
tion concerning damping ratio and natural frequency may be 
obtained after the plot is complete by examining individual 
points, as described in the previous paragraph. Figures A.1 
and A.2 are examples of DACSAP root locus plots. 

It should be noted that if one or more lines of heading 
are specified, the flot will be rotated on its side so that 
the long edge of the plot is aligned with the long edge of 
the screen. With the plot oriented in such a way, a hard- 
copy print of the screen will tretter fill the page. eae eC 
user intends to run a multitude of plots, conducting anal- 
ysis cn the screen, it 1s suggested that the plots be 
produced without headings until a hardcopy is desired. This 
will reduce the plotting time and produce a larger, upright 
plot to analyze on the screen. 

Subsequent selection of Option 1 from the Root Locus 
Options menu will result in the following menu of change 
options: 

| ROOT-LOCUS CHANGE OPTIONS =| 


a I EE ee el 


| RANSFER FUNCTION(S) 
2 AIN PARAMETERS 
3 CHANGE PLCTTING PARAMETERS 
4 CHANGE PLOT HEADING 
5 NO CHANGES (READY TO PLOT) 
6 | RETURN TO DACSAP MENU 
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Figure A.1 Sample w'-plane Root Locus Plot. 
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Sample z-plane Root Locus Plot. 


Figure A.2 
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From this menu the user may selectively change any or all of 


the parameters required to produce another root locus plot. 
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De. FREQUENCY RESPONSE ANALYSIS 


oe SD OD De ee es ee ee es ee es ee es ee ee ee ee eee ee ee ee eee ee ee eee ee 


| OPTICN NO. ] OPTION | 
| 1 | INPUT TRANSFER FUNCTION (S) | 
2 | ROOT LOGUS ANALYSIS | 
? 5 BODE ANALYSIS (OPEN LOOP) i < 
> 4 BODE ANALYSIS (CLOSED LOOP) } < 
> 5 NYQUIST ANALYSIS 1} < 
> | 6 | NICHOLS ANALYSIS | < 
7 | TIME RESPONSE | 
8 CHANGE BLOCK DIAGRAM | 
5 SAVE THAIS PROBLEM 
10 START A NEW PROBLEM 
17 HELP 
12 BxLlDACSAPL | 
DACSAP offers four different freguency domain analysis 
mools. They are open and closed loop Bode plots, polar 


plots (for Nyquist analysis) and log magnitude-phase plots 
(for Nichols analysis). ODtLONS Semeur Ough Gaot the DACSAP 
Options menu should be used to select the desired analysis 
technique. 

The Frequency Response Options menu, below, is presented 


when any of the four frequency domain methods are chosen, 


a i co SE cS ei Sc Sc ce ec ee ce em ee ee 


| OPTION NO. | OPTION 
1 | GRAPHICAL ANALYSIS 
2 TABULATED DATA 
3 RETURN TO DACSAP MENU | 
U EXIT DACSAP | 
5 | HELP 


For the three open loop analysis methods (Bode, Nyguist 


and Nichols) the following information will be requested: 


mae OUnNCY RESPONSE PARAMETERS: 
WHAT IS YOUR OPEN-LOOP GAIN? 
GAIN = 


FREQUENCY RANGE: 
MIN FREQUENCY (RAD/SEC) 
MAX FREQUENCY = (RAD/SEC) 
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PLOT HEADING 
HOW MANY LINES OF HEADING WOULD YOU LIKE? (4 5ieaee 
THON OF 22 Cece tay 
The frompts for closed loorp Bode analysis are tne same 
except that instead of OPEN LOOF GAIN the user is asked for 
the following: 


WHAT IS THE FORWARD PATH GAIN? 
WHAT IS THE FEEDBACK PATH GAIN? 
POSITIV=2 OR NEGATIVE FEEDBACK? (P OR N) 


The more detailed information requested for closed loop Bode 
plots is required so that the closed loop transfer function 
may be calculated. — 

Tf a discrete time system is being analyzed, the 


following question will also be asked: 


FICTITIOUS (V) OR REAL {W) FREQUENCIES? (VV Ck 


Since w- and w't-transformations result in frequency 
"warping", a frequency transformation must take place if the 
response is to be plotted as a function of £real frequedag 
DACSAP will perform the transformation if W (omega) 1s 
selected or plot the response as a function of fictitious 
frequency if V (nu) is chosen. If Wis selected the user 
will be prompted for the sample period which is required to 
perform the frequency transformation. 

Once all of the parameters have been input, the 
frequency response will be plotted using the type o£f plot 
selected from the DACSAP Option menu. Figures A.3, A.4 and 
A.5 are examples of DACSAP frequency response plots. 
Frequency response plots are automatically scaled, so the 
user is not prompted for plot dimensions. 

When these plots are complete, information concerning 
phase and gain margins and crossover freguencies is 


presented as follows: (sample values) 
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Figure A.4 Sample Polar Plot. 
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-375-850.0-325.0-300.0-275.0-250.0-225.0-200.0-175.0-150.0-125.0-100.0 -75.0 -50.0 -25.0 
PHASE (deg) 


Figure A.5 Sample Log Magnitude-Phase Plot. 
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] DEGREE |] NUMERATOR | DSNOMINATOR | 

| 4 ] ] 1.0000000 | 
3 j { 9.0400000 | 
2 25-0000 000 QO. 3760009 

1 25.7499 847 0. 20 7/9999 

| 0 | O2 7999999") 0.5/599 20a 

PHASE CROSSOVER FREQUENCY = 0.5299910E+00 RAD/ZSuES 

GAIN MARGIN = —24 858267 DE 

GAIN CROSSOVER FREQUENCY = 0.4135011E+01 RAD/SEC 

PHASE MARGIN = 51.0 15 50D rare rs 


As with root locus plots, folar and log magnitude-phase 
plots are rotated on their sides if one or more lines of 
heading are specified for the plot. 

Subseguent selection of Oftion 1 from the Frequency 
Response Options menu will result in the following menu of 


Change options: 


| OPTION NO. | ‘ OPTION | 
1 CHANGE TRANSFER SNe a Cone | 

Z CHANGE FREQUENCY RESPONS& PARAMETERS | 

ss | CHANGE PLOT HEADING | 

a NO CHANGES (READY TO PLOT) | 

| > | RETURN TO DACSAP MENU j 


=P cum cp ee ae SD we cs ee ee ees es ee es eee ee ee ee ee eee ee es ee ee eee eee ee ee ae ee eee ee 


From this menu the user may selectively change any or all of 
the parameters reguired to produce another frequency 


response plot. 
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Eo TIME RESPONSE ANALYSIS 


ee ee ee es ee ee ee ee ee ee ee ee ee ee oe 


Poel ON NOe| OP Eon | 
1 | INPUT TRANSFER FUNCTION (5S) | 
Ze | ROOT LOCUS ANALYSIS | 
: BODE ANALYSIS (OPEN LOOP) | 
4 BODE ANALYSIS {CLOSED LOOP) | 
3s NOU 2S TAN AInY ses | 
6 Niet OVsS wa vALY SES | 
> y/ Poe ek SoPONS Ss i < 
8 CHANGE ELOCK DIAGRAM { 
2 SAVE elas PROBLEM | 
10 START A NEW PROBLEM { 
11 | HELP | 
12 EXIT DACSAP { 


A time history of the total system's response to a user 
specified input may be obtained by selecting Option 7 from 
the DACSAP Option menu shown above. This selection will 
result in the presentation of the following menu of time 


response options: 


a ce ce ee ce ee es ee ee ce ec ee ee ce ee es es ee es ee ee ee ee ee ee eee 


| OPTION NO.] OPTION l 
1 GRAPHICAL ANALYSIS ' 

2 TABULATED DATA 

3 RETURN TO DACSAP MENU 

y EXIT DACSAP | 
5 HELP 


Selecting Options 1 or 2 will result in the user being 


prompted for the following parameters: 


Pont RESPONSE PARAMETERS: 
1 52 BP 
2 pipe Ue San 
3 = RAMP 
wee ser TYPE OF DPNPUT TO YOUR SYSTEM (1,2 OR 3) 
WHAT IS THE AMPLITUDE OF YOUR INPUT? 
Weat iS THE FORWARD PATH GAIN? 
MEAd Do THE FEBDBACK PATH GALN? 
POSITIVE OR NEGATIVE FEEDBACK? U2 eB 


WHAT IS YOUR SAMFLE PERIOD? T = {SEC) 
(discrete systems) 
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FOR HOW MANY SECONDS DO YOU WISH TO SEE THE RESPGiS 


If Option 1 was selected from the Time Response Options 
menu, the user will also be frompted for a plot heading. 
Time response plots are automatically scaled, so the user is 
not prompted for plot dimensions. If Option 2 was chosen, 
the output may be selected to go to the screen or the 
PEIiDteLrs Figures As«6 and A.7 are examples of DACSAP time 
response plots. 

Subsequent selection of Option 1 from the Time Response 


Options menu will result in the following menu of change 


options: 
| TIME RESPONSE CHANGE OPTIONS | 
| OPTION NO. | QE LON { 
1 | CHANGE TRANSFER FUNCTION (S | 
2 CHANGE TIME RESPONSE PARAMETERS 
S| CHANGE PLOT HEADING 
4 | NO CHANGES (READY TO PLOT) | 
| 5 |] RETURN TO DACSAP U | 


From this menu the user may selectively change any or all of 


the parameters required to produce another frequency 
response plot. 
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F. CHANGING THE BLOCK DIAGRAM 


inieOoLdcr tom change of add to the transfer functions in 
the block diagram, Option 8 of DACSAP's option menu should 


be selected. 


—_ <= =P ee ee ee et eee ee ee ee A ee ees ee ee ee ee eee eee eee eee eee eee eee ee eee 
a a ee ee ee ee ee ee eee eee eee eee eee eee ee ee eee eee 


OPEN LOOP) 
: CLOSED LOOP) 


E 
CHANGE ELOCK DIAGRAM 
SAVE THIS PROBLEM 
een A NEW PROBLEM 


EXIT DACSAP 


Entry to this change routine can also be obtained by 
selecting Option 1 of the root locus, freguency response or 
time response parameter change menus. 

Transfer function change oftions are presented in the 
following menu: 


{ cio foneruNnel LON s(T .F =) OC CHANGE OPTIONS { 


] OPTION NO. | OPTION | 
i | CHANGE A T.F. IN THe CURRENT LOOP l 
Z ADD A UNE WeBECCK TG THE CURRENT LOOP jf 
3 {| CHANGE T.F. IN AN INNER LOOP { 
Ly EXPAND TO AN OUTER LOOP 
5 NO CHANGES 
6 BEEP 


iaewicy current, Loop 


When Option 1 of the Transfer Function Change menu 


is selected, the user will be asked which block is to be 


changed. Then the following menu of block change options 
wili be presented: 


BLOCK CHANGE OPTIONS 
BLOCK NUMBER 1 


| OPTION NO. | OPTION | 
1 CHANGE CURRENT NUMERATOR 
2 CREATE A NEW NUMERATOR 
3 CHANGE CURRENT DENOMINATOR 

| 4 CREATE A NEW DENOMINATOR | 
5 CREATE NEW NUMERATOR & DENOMINATOR | 
6 MOVE BLOCK 10 THE FEEDBACK PATH | 
7 NC MORE CHANGES TO THIS BLOG: | 


If Options 1o0r 3 are selected, asummary of the 
roots or coefficients of the selected polynomial will be 


presented in one of the forms shown below: 


= SD wb eS Se SS a Se SS ee SE Se SS SP ee ce ee SE eG Se 


| LOOP NUMEER 1 | 
BLOCK NUMBER 1 | 
| LTE NO. | NUMERATOR ROOTS | 
1 | -0.0300000 0.0000000 J 
2 CONSTANT = 25. 9000000 
ANY CHANGES TO THESE PARAMETERS? 


Se cee 


| LOOP NUMBER 3 | 
BLOCK NUMBER 2 | 


mm ee cc eee wc ee ee ee ee 


| TEM No. | NUMERATOR COEFFICIENTS | 
1 1.0000000 S$ ** 2 ] 
Z =0 236000000 53° ae | 
8 0. 1600000 S>== | 
ANY CHANGES TO THESE PARAMETERS? 


By answering yes (Y) to the question following the summary, 
the user may change any of the parameters displayed. A 
negative response will return the user to the Block Change 
Options menu. This method may be used if the user simply 


wishes to review the contents of a block. 
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If Options 2 or 4 are selected from the Block Change 
Options menu, the user will be prompted for all new parame- 
ters for the selected polynomial. The form of the polyno- 
mial (factored or coefficient) will remain the same as it 
was before. However, all other parameters may be changed 
including the order of the polynomial. The user will be 
prompted for the new parameters as he was in the transfer 
function input routine. 

If the user wishes to chanye the entire character of 
a block, Option 5 of the Block Change menu should be 
selected. This option essentially erases the current paran- 
eters associated with the selected block and prompts’ the 
user for all new parameters. The user is free to input the 
transfer function inany form and entry may be from _ the 
console or from a data file. The prompts are exactly as 
described in the Transfer Function Input section of this 
manual. 

Option 6 of the Block Change Options menu simply 
allows the user to redesignate a block as being in the feed- 
back path if it is currently in the forward path or vice 


versa. Option 6 will read 
| 6 | MOVE BLOCK TO THE FORWARD PATH | 


for blocks which are currently in the feedback path. 
When no more changes to the currently selected block 
are required, Option 7 should be selected. The user is then 


returned to the Transfer Functicn Change menu. 


V7 


— a ee ee eee eee eee eee ee eee ee PP See ee eee ee ee ee eee ee eee ee ee ee ee ee ee eee ee ee eee ee 
= <> ee oe ee ee oe es ee ee ee es ee ee es es ee es es es eee ee ee ees ee es es eee es es ee ee eee ee 


. IN THE CURRENT LOOP l 
O THE CURRENT LOOP | 

INNER LOOP | 
R LOOP | 
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By selecting Option 2 from the Transfer Function 
Change menu, above, the user is simply sent to the transfer 
function input routine to enter the parameters for the new 
block. The prompts are exactly as described in the Transfer 
Function Input section of this manuai. The transfer func- 
tion may be of any form, it may be placed in the forward or 
feedback path and may be entered from the console or froma 
data file. 


If the user is currently working on an outer loop of 
a multiloop system where the inner loop(s) were designed and 
saved using DACSAP, then the inner loop(s) may be returned 
to and changed. Later it will be shown how to have these 
changes reflected in the outer loop(s). 


] TRANSFER FUNCTION (T.F.) CHANGE OPTIONS | 


| OPTION NO.| OPTION 
1 CHANGE A T.F. IN THE CURRENT LOOP 
2 ADD A NEW BLCCK TO THE CURRENT LOOP | 
> 3 CHANGE T.F. IN AN INNER LOOP [ea 
| 4 EXPAND TO AN OUTER LOOP 
5 | NO CHANGES ' 
6 | HELP 


When Option 3 of the Transfer Function Change 


Options menu is selected the user is prompted as follows: 


70 


Pees OU AREM@GOLNG 70 WANT TO COME BACK TO THIS @600P 
oon, LOU OUDE SAVE IT: 


PO OUSWISH TD) PSAVE THE TRANSFER FUNCTION(S) IN 
THIS LOOP? (Y OR QW) 

As stated, the current loop must be saved if it is to be 

returned to after changes are made to tne inner loop. ff it 

is not saved, its transfer functions will have to be 

re-entered from the console. The following subsection will 

cover how to get back to the outer loop(s). 

When the user selects to save the current ioop he 
will be prompted for a filename under which it is to be 
saved. After the outer loop is saved, the user will be 
prompted for the filename of the inner loop, that loop will 
ke loaded and the following guestion will appear at the 
terminal: 

feeiedo MAKE ANY CHANGES TO YOUR TRANSFER FUNCTION (S) ? 
(Y OR N) 

If a negative response 1S input, DACSAP's main 
option menu will be presented and an analysis technigue for 
the inner loop may be selected. An affirmative response 
will result in the Transfer Function Change menu being 
presented and the user is free to change or add to the inner 
loop's transfer functions. It should be noted that if 
nothing but the forward or feedback path gains are to be 
changed, Option 4 of the Transfer Function Change menu 
should be selected to re-expand to the outer loop. The 
Opportunity to change these gains will be given at that 


time, as explained in the following subsection. 


4. Expanding to an Outer Loop 


If the system currently being worked on is an inner 
loop of a muliti-loop system, the user may expand the svsten 


to an outer loop, letting DACSAP calculate the inner loop 


1 


closed loop transfer function and place that transfer func- 
tion anywhere in the cuter loop. This operation can be made 
to occur by selecting Option 4 from the Transfer Function 


Change menu. 


= a ee ee ee eee ee ee ee ee ee ee eee ee ee ee ee ee ee ee ee ee ee eee ee 
= oe SP ee ee ee Oe eee ee eee eee eee ae AE Eee EE ae aE ewan 4A aE aE eRe 4A Ea ae 2 2p SER SP ae ee eS eee ae ae 


1 | CHANGE A T.F. IN THE CURRENT LOOP | 
Z ADD A NEW BLCCK TO THE CURKENT LOOP | 
3 CHANGE T.r. IN AN INNER LOOP | 
> 4 EXPAND TO AN OUTER LOOP < 
> | NO CHANGES 
6 HEE 


Before expansion to the outer loop occurs, the user 
is given the opportunity to save the inner loop to a data 
Hae r. This should be done if the inner loop is ever to be 
returned to for further analysis or changes. Prompting will 
be as follows: 

IF YOU THINK YOU MIGHT WANT TO GO BACK AND MAKE 
CHANGES TO THE CURRENT LOOP, YOU SHOULD SAVE IT. 


DO YOU WISH TO SAVE THE TRANSFER FUNCTION(S) IN 
THIS LOOP? (Y OR N) 


WHAT NAME DO YOU WISH TO GIVE TO YOUR DATA FILE? 
(8 CHARACTERS MAX) 


| CLOSED LOOP PARAMETERS | 


{1 ITEM NO. | PARAMETER 1 VALUE | 
1 FORWARD PATH GAIN 1.000] 
Z FEEDBACK PATH GAIN 1.57 
3 POS. OR NEG. FEEDBACK N { 


SSS SS SSS SE ee eee ee ee ee ee ee ee ee ee ee eee ee 


THE CLOSED LOOP TRANSFER FUNCTION WILL ALSO BE 
CALCULATED AND SAVED. 


ANY CHANGES TO THESE PARAMETERS? 


In the case shown above, the closed loop parameters had 
already been defined during the analysis of the inner loop. 
So, the user 1s Simply given a last opportunity to change 
them. If they had not already been defined, the user would 
be prompted for them at this time. 
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MeeeiismpOllcGmett 2S time to define the outer loop. 


The following prompt will begin the process: 


DO YOU WISH TO INPUT THE OUTER LOOP TRANSFER FUNCTION (S) 
FROM A DATA FILE OR FROM THE CONSOLE? (D OB C) 


Entry from the console will be discussed first. 


a. Entering the Outer Loop From the Console 


The first time that the system is expanded, the 
outer loop must be entered from the console. The user will 
be prompted as follows: 

INTO WHICH BLOCK OF THE NEW LOOP DO YOU WISH TO PLACE 
fi CLOSED LOOP TRANSFER FUNCTION OF THE INNER LOOP? 
Momititis BLOCK IN DHE FORWARD OR FEEDBACK PATH? {F OR 2B) 


fie ClLOSED LOOP TRANSFER FUNCTION FROM THE INNER 
Peve SAS BEEN PLACED IN BLOCK NO. 2 OF THIS LOOP. 


HOW MANY BLOCKS (INCLUDING ELOCK NO. 2) ARE IN 


THIS LOOP 
In this case, the answer to the second question was 2. The 
system may have looked like the one below. So, the closed 


loop transfer function for the inner loop (enclosed by the 
dashed line) becomes block number 2 of the outer loop. FOr 
this systen, the response to the last question would be 4. 

The user 1S now prompted for the transfer func- 
tion parameters for the remaining blocks in the outer loop 
fELocks 1, 3 and 4 in the example). These prompts are 
exactly like those presented in the Transfer Function Input 
section of this manual. 

It 1s very important to note that the block 
Number containing the inner Jloop's closed loop transfer 
Eunction (block number 2 in the example) is always retained 
as a permanent parameter associated with an outer loop. 
This becomes important whenever an inner loop is changed and 
the outer loop is subsegquently returned to, as discussed in 


the next subsection. 


8 1 


= 
BLOCK 3 





Figure A.8 Example Systen. 


b. Entering the Outer Loop From a Data File 


This option should te used whenever the user has 
returned to an inner loop to make changes and is now ready 
to re-expand to the outer loop. The outer loop must have 
been saved before the inner lcop was loaded in order for 
this oftion to be used. 

The first time that the system was expanded to 
the outer loop, the user specified a block in the outer loop 
into which the inner would be placed. That block number 
became a permanent parameter associated with the outer loop. 
For this reason, when changes have been made to the inner 
loop, they will automatically be reflected by a change to 
the correct block in the outer loop. 

After the user has chosen to expand to the outer 
loop (using Option 4 from the Transfer Function Change menu) 
and defined the closed loop parameters for the inner loop, 
the following prompts will be presented: 


DO YOU WISH TO INPUT THE OUTER LOOP TRANSFER FUNCTION (S) 
FROM A DATA FILE OR FROM THE CONSOLE 2 Dee rae. 
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Mintel S “HEeNAMes, OF YOUR DATA ELLE? (FILE NANE ONLY) 


Pome LOokD POOP RIRANSTER FUNCTION FRO THE INNER 
LOOP HAS BEEN PLACED IN BLOCK NO. 2 OF THIS LOOP. 


Naturally, one would respond with aD (data 
file) to the first question then provide the filename for 
the outer loop. The program then loads the outer ioop 
transfer functions, but replaces the transfer function in 
the expansion block (rlock number 2 in the example) with the 
updated closed loop transfer function for the inner loop. A 
new open loop transfer function for the outer loop is also 
calculated using the new closed loop transfer function held 
in the expansion block. The user is now asked 

WANT TO MAKE ANY CHANGES TO YOUR TRANSFER FUNCTION (S) ? 
(Y OR N) 
An affirmative response will send the user to the Transfer 
Function Change menu for the outer Loop. A negative 
response sends the use to the LACSAP Options menu where he 
May select an analysis technicue and observe how the changes 
to the inner loop have affected the total system's 


performance. 


5- Manipulating Large Multiloop Systems 


Much of the power of DACSAP is due to its ability to 
Manipulate large multiloop systems through the use of 
Options 3 and 4 of the Transfer Function Change menu. The 
key to making these features work is to insure that the 
system is input starting with the innermost loop’ then 
expanding one loop ata time, Saving each loop along the 
way. If this is done then even the innermost loop may be 
returned to and changed having those changes automatically 
incorporated in the outer loop (via Option 4 of the Transfer 


Function Change menu). It should be noted that the system 


es 


must be re-expanded one loop at a time when changes have 
teen made to an inner loop. This is so each successive loop 
may be updated with the change to the inner loop. 

The process 1s made particularly easy if the user 
Simply wishes to make changes to the gains of an inner loop, 
which is often the case. After the inner loop is loaded, 
Simply select Option 4 of the Transfer Function Change menu 
to re-expand to the outer loop. The gain parameters for the 
inner loop will then be presented and may be changed before 
the outer loop is loaded and updated. 

What if an outer loop has 2 (or more) ianer loops in 


series, as shown in Figure A.9? This situation and others 





Figure A.9 Inner loops in Series. 


like it can be handled by designing the inner loops as 
Subsystems and saving them to data files. Tneo when 
designing the outer loop, these subsystems may be loaded 
into their appropriate blocks. If these subsystems must 


later be changed, then the outer loop should be saved and 
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the desired inner loop loaded as a new problem (Option 10 of 
the DACSAP Options menu). The inner loop should then be 
Changed, re-saved andthe outer loop loaded aS a new 
problen. Then the old subsystem may be replaced by the new 
one through Option 5 of the Block Change menu. 

The user may end a terminal session and, later, pick 
up where he left off if he saves the system before exiting 
the progran. If the user has designed and saved an inner 
loop and, during a new terminal session, wants to begin work 
on the next loop (Loop 2, for example), he may do so. 
Recall that, during the transfer function input sequence, 
the user is asked for the number of the current loop. In 
this case, the user would respond with a2 and would be 
asked for the block number in Loop 2 which contains’ the 
closed loop transfer function fcr Loop 1. When loading the 
designated block, the user will be prompted for the filename 
under which the inner loop (Loof 1) is saved. Changes to 
the inner loop may then be made using Options 3 and 4 of the 
Transfer Function Change menu. 

One can see that loops and transfer functions may be 
manipulated to form systems of fractically any shape. Very 
complex systems may be designed through the imaginative use 
of the subsystem block load and loop expansion features of 
DACSAP. 


85 


G. SAVING THE SYSTEM TO A DATA FILE 


| OPTION NO. | OPTION | 
1 INPUT TRANSFER FUNCTION (S) | 
Z ROOT LOCUS ANALYSIS | 
5 BODE ANALYSIS fares Looe | 
4 BODE ANMEI STS (CLOSED LOer 
| S | NYOUTS tea Ney Sls { 
6 NICHOLS Bane | 
7 | TIME RESPONS | 
8 | CHANGE ELOCK DIAGRAM { 
> 9 SAVE THIS PROBLEM | < 
a START A NEW PROBLEM 
eZ EXIT DACSAP 
Option 9 from the DACSAP Option menu, above, allows the 
user to store the system's parameters to a data file. This 


may be done so that the system can be loaded during a 
subsequent terminal session for further analysis or to 
continue the design frocess. The user will be presented 
with the following two questions: 

WHAT NAME DO YOU WISH TO GIVE TO YOUR DATA FILE? 

(8 CHARACTERS MAX) 


DC YOU WISH TO HAVE THE CLOSED LOOP TRANSFER 
FUNCTION CALCULATED AND SAVED? (Y OR N) 


The response to the second guestion iS important. Recall 
that when inputting transfer functions, any block may be 
loaded from a data _ file. The transfer function that is 


loaded at that time depends on which one is saved now. 

If the system being saved is a feedback system, then the 
closed loop transfer function should be calculated and 
saved. Then it can be loaded as a subsystem into a block of 
a larger systen. 

If the system being saved is a compensator, filter or 
any other system consisting of a Single component or compo- 
nents in cascade, then the closed loop transfer function 


should not be calculated. In this case, the cascade 
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transfer function will be saved and may be loaded into a 
block of another system. Note that the system will be saved 


on the user's 'A‘t disk as 
SEP EPeENAME. | wDAT A Al 


where 'FILENAME' is the name specified by the user. 


H. THE HELP ROUTINE 


A HELP routine exists in DACSAP which contains useful 
information pertaining to all of the program's najor 
routines. These HELP files contain a list of required input 
data, brief operating instructions anda brief description 
of options for each routine. HELP may be selected fron 
DACSAP's main option menu, the transfer function input and 
change henus, and each of the analysis routines' option 
menus. If HELP is selected from DACSAP'sS main option nenu, 
any of the HELP files may be examined or the entire HEL? 
routine may be printed out. If HELP is selected from any 
other option menu, informaticn pertaining only to that 
routine is presented, then the user is returned to the menu 
from which HELP was called. 
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ot DAC Steen anoPokR PUNGCTIONSINPUT 2OCTINE *** 


PURPOSE: 

THIS ROUTINE INPUTS THE TRANSFER FUNCTION DATA an 
meee ALi OF DACSAPS ANALYSIS ROUTINES. Te Lip ieewALA NUST 
BE AN OPEN-LOOP TRANSFER FUNCTION WHICH WILL BE INPUT AS A 
RATIO OF POLYNOMIALS DESCRIBED BY EITHER POLYNOMIAL 
COEFFICIENTS OR POLYNOMIAL ROOTS (FACTORED FORM). 


REQUIRED DATA: 


1. TRANSFER FUNC CTION (S) WHOSE EOLYNOMIALS ARE IN FACTORED 
OR COEFFICIENT FORM 

2. COEFFICIENTS MUST BE REAL NUMBERS ENTERED IN DESCENDING 
POWERS OF S, Z, WOR ¥! 

3. IN FACTORED FORM, THE COMPLEX ROOTS OF THE POLYNOMIAL 
MUST BE ENTERED ALONG WITH A REAL GAIN CONSTANT 
(MULTIPLICATION FACTOR) . 

PROCEDURE: 


fee OhrriCIENT FORM; 
GIVEN THE POLYNOMIAL 
Zesty tai ODOn 5) to 
THE SCOETOIVCEENTS WOULD BE ENTERED EN THE FOLLOWING ORDER 


2 
e050 
8 


A PROMPT WILL APPEAR BEFORE EACH ENTRY. 
2. FACTORED FORM; 
THE SAME POLYNOMIAL MAY BE REPRESENTED AS 
PC ae so ome omen teuq th — 1.0443) 
AGAIN, THE USER WILL BE PRCMPTED. INPUT DATA WILL BE 


CONSTANT = 2 
ROOT 1: REAL PART = -1.414 
IMAGINARY FART = -1.414 
ROOT 2 REAL PART = -1.414 
IMAGINARY PART = 1.414 
Semel FOLLOWING RESTRICTIONS APPLY TO THE TRANSFER 
MUNeTIONS; 
Poe nO Um bLOCKS NAY BE ENTERED 2NTO THE CONTROL LOOP. 
Peeeonen oLOCK MAY BE UP TO 20TH ORDER. 
C. ane oe Cretns LOTAL > tslEteeuSteNOL EXCEED 100TH 
De. THE NUMERATOR AND DENOMINATOR OF ANY GIVEN BLOCK HU 
Poot lease FORM (CORFFICIENT OR cee he 
Deeley sel se LOCK IS ENTERED GS ITIieiS ASSUMED [TO BE IN 
THE FORWARD PATH FOR CLOSED-LOOP CALCULATIONS. 
Poon y FEEDBACK ES ASSUNED TF NG BLOCKS ARE ASSIGNED 
TO THE FEEDBACK PATH. 
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OPTIONS 7 COMMS: 


1. ANALYSIS MAY 
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~et DACSAPEIRALNor ER FUNCTION CHANGE ROUTINE *** 


PURPOSE: 

Peso OULINE PROVIDES THE USER WiTH THE ASLTLITY TO eXAAINe 
AND / OR CHANGE TRANSFER FUNCTIONS IN THE CONTROL LOOP 
CURRENTLY BEING DESIGNED OR ANALYZED. THE ROUTINE ALSO 
Bee ws nt USER TO EXPAND TO OUTER LOOPS OR RETURN TO AND 
CHANGE INNER LOOPS. 


REQUIRED DATA: 


ae POLENOMTAL (S) WHICH ARE TO BE CHANGED OR ADDED WILL BE 
INPUT IN THE SAME MANNER AS DESCRIBED IN THE TRANSFER 
FUNCTION INPUT SECTION. 

2. IF EXPANDING TO OUTER LOOPS, THE INNER LOOP DATA MUST BE 
SAVED IF THE USER EVER INTENDS ON RETURNING TO AND 
CHANGING ANYTHING IN PHE INNER LOOP. THE USER WILL BE 
PROMPTED FOR A FILE NAME DURING THIS PROCEDURE. 

3. FORWARD AND FEEDBACK PATH GAINS AND THE TYPE OF FEEDBACK 
ARE REQUIRED WHEN LOOPS ARE SAVED. THIS DATA IS NEEDED 
IN ORDER TO CALCULATE THE LCOPS CLOSED LOOP TRANSFER 


PROCEDURE: 


ie BERGE ieee COPPA EON FROMTHE SERANST ER FUNCTION 
A I 
Zee LHE USER Wie, PROWET SD SV HTLE sia RENG SGHANGES TO ANY 


cele THe, FORM OF AN EXISTING BLOCK IS TO BE CHANGED, THE 
Sin MisteseLECT THE OPTION TO CREAT® SOTH A NEW 
UMERATOR AND DENOMINATOR FOR THAT BLOCK. 
Mea eN KRE-ENTERING INNER LOOPS, SAVE THE CURRENT LOOP AND 
ROVIDE THE FILE NAME O HE LOOP TO BE RE-EXAMINED. 
XPANDING BACK OUT TO T OUTERSUOSieLO@P UUST BE DONE 
Nimet@Or Ast A TIME. TH WILL ALLOW THE PROGRAM TO 
DHOMATICATLY UPDATE THE OUTER LOOPS IWITH THE CHANGES 
ADE IN THE INNER LOOP. ; 
Seo AVE BACH LOOP DURING THE EXPANSION PROCESS SO THAT THE 
POerewtLbs CONTAINS THE MOST UP=-TO=DATE INFORMATION. 


OEE ONS / COMMENTS; 


L 
if 
HE 
LS 
E 


ret Otity =tic 


ieee AEOLING TRANSFER FUNCTIONS MAY BE CHANGED OR MOVED. 

Ze Bee heNSE Es FUNCTION(S) MAY Be ADDED TO THE CONTROL 

3. THE CONTROL SYSTEM MAY BE EXPANDED TO MULTIPLE LOOPS. | 
iti NeR LOOP CLOSED LOOP TRANSFER FONCTION IS AuUTC- 
woop Heer ert O 2 USER pEE PED BLOCK EN THE OUTER 

4. INNER LOOPS WHICH HAVE BEEN SAVED MAY BE LOADED AND 
mainvGob. ONCE CHANGES HAVE BEEN MADE, THE OUTER LOOPS 
CAN AUTOMATICALLY BE UPDATED iF EXPANSION IS DONE ONE 


LOOP AT A TIME. 
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***X DACSAP ROOT LOCUSRROUIING == 


FURPGsSr. 

THIS ROUTINE CALCULATES THE SYSTEMS CLOSED LOOP) CE ARAGa 

ISTIC EQUATION FROM A PREVIOUSLY ENTERED OPEN LOOP Tk 

FUNCTION. THE ROOTS OF THE EQUATION ARE DETERMINED Olea 

USER DEFINED RANGE OF OPEN LOOP GAIN VALUES AND MAY BE 

PLOTTED AND/OR PRESENTED IN TAEFULAR FORM. 

REQUIRED DATA: 

1. A PREVIOUSLY ENTERED OPEN LOOP TRANSFER rUNCTION. 

2. A RANGE OF OPEN LCOP GAIN VALUES. THESE GAINS WAY BE 
POSITIVE OR NEGATIVE, BUT NCT BOTH. 

3. TYPE FEEDBACK AND SAMPLE PERIOD {DISCRETE Sis 7fi ae 

4. PLOTTING DIMENSIONS. RECTANGULAR DIMENSIONS 
ee XHAX-AMIN = EO Se aene SHOULD BE USEDL IF ANGI Seer 
MAGNITUDES WILL BE MEASURED DIRECTLY FROU THESE 
IF UNSURE OF THE PLOTTING DIMENSIONS, SELECT TABULAR 
DATA FIRST. POLE AND ZERO LOCATIONS WILL BE PROVIDED ese 
THAT PLOTTING DIMENSIONS MAY BE DETERMINED. 

PROGECU RE: 

1. SELECT PLOTTED OR TABULAR DATA OPTION. 

2-« INPUT GAIN AND PLOTTING PARAMETERS WHEN PROMPTED- 

3. AFTER THE ROOT-LOCUS HAS BEEN PLOTIED, THE GAIN, Nati 
FREQUENCY, AND DAMPING RATIC FOR SPECIFIC POINTS ON THE 
LOCUS MAY BE REQUESTED. THEN, THE ROOTS ASSOCIATS Daa 
ANY GAIN VALUE MAY BE HIGHLIGHTED ON THE PLOT. 

+. TABULAR DATA MAY EE OBTAINED AFTER PLOTTING. 

OPTIONS 7 VCOnMEN TS: 


1. TEMPLATES FOR CONSTANT DAMPING RATIO AND NATURAL 
FREQUENCY ARE PROVIDED FOR DISCRETE DOMAIN 2lOle 

Za GAIN, DAMPING RATIO, AND NATURAL FREQUENCY MAY BE 
REQUESTED FOR ANY POINT ON THE ROOT-LOCUS. 


3. A DESIGN POINT MAY BE HIGHLIGHTED ON THE PLOT BY 
SPECIFYING ITS GAIN VALUE. 
4. O-4 LINES OF HEADING MAY BE SPECIFIED FOR EACH PLOTS 


a2 


<i veowaeentOurNCY RESPONSE ROUTINE =+* 


EURPO SE : 
But o ROUTINE CALCULATES AND PLCTS THE eae Feo oLONSss OT 
Oe EN OR CLOSED SDOGOr SIYSTEN JO A USER S2ECIFIED INPUT 
SIGNAL AND FREQUENCY RANGE. OPEN LOOP RESPONSE MAY BE 
mreebi ZED USING BODE PLOTS, PCLAR PLOTS et ee ANALYSIS), 
On PHASE-SAGNITUDE FLOTS Ve ANALYSIS). €HOSED LOGE 
RESP ONS® MAY BE ANALYZED USING BODE PLOTS. 
REQUIRED DATA: 
eee PREVIOUSLY ENTERED OPEN LCOP TRANSFER FUNCTION. 
2- OPEN LOOP GAIN (BROKEN INTO FORWARD AND FEEDBACK PATH 
oA NS@2eR Cles Eo eECOP RESPONSE) 
Pepa ANGE OF FREQUENCIES. 
4. SAMPLE PERIOD (DISCRETE SYSTEMS). 
Semeexy PE OF FEEDBACK (DISCRETZ SYSTEMS) 
PROCEDURE: 
t. SELECT PLOTTED OR TABULATED DATA OPTION. 
2- INPUT THE REQUIRED DATA (SEE ABOVE) WHEN PROMPTED. 
feet, PLOTS ARE AUTOMATICALLY SCALED. FOR NYQUIST PLOTS, 
ira LARGE FREQUENCY RANGE I5 SELECTED, THE AUTOMATIC 
SCALING CAN RESULT IN LOW RESOLUTION NEAR THE ORIGIN. 
PNCKREASDING THE MIN FREQ MAY RESULT IN A BETTER PLOT. 
4. TABULAR DATA MAY BE OBTAINED AFTER PLOTTING. 
OPTIONS / COMMENTS: 
1. OPEN LOOP PHASE AND GAIN MARGINS AND CROSS-OVER 
PREQUENCIES ARE FROVIDED AFIER EACH PLOT. 
peer on THE DISCRETE DOMAINS, FREQUENCY WARPING OCCURS 
DURING THE CALCULATIONS. THE UGSER SPECIFIES WHETHER 
PiomenrolONoe £5 10 BE PLOTTED AGAINST THE FICTITIOUS 
ene ont FREQUENCIES OR REAL FREQUENCTES. 
comes LINES OF HEADING MAY BE SPECIFIED FOR SEACH PLOT. 
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**¥* DACSAP TINE RES RCTS ROU Tr ae 


PURPOSES. 


THIS ROUTINE CALCULATES AND PLOTS THE TRANSIENT RiSeroleae 
THE CLOSED LOOP SYSTEM 70 AUCH PoP eCiit i hee 
ROUTINE CALCULATES THE CLOSED [LOOP TRANSFER FUNCTION Ves 
OPEN LOOP DATA AND INPUTS TO Sapo eroOUr Tee 


REQULRED Data: 


1. A PREVIOUSLY ENTERED OPEN LCOP TRANSFER FUNGI Tee 
- TYPE AND MAGNITUDE OF THE INPUT SIGNAL. 
3. FORWARD AND FEEDBACK PATH GAZTNS AND THE TYPEReT 
FEEDBACK. 
G4. SAMPLE PERIOD see Chae a oe 
5. THE NUMBER OF SECCNDS OVER WHiCH THE RESPONS?: Dome 
OBSERVED. 
PROCEDURE: 
1. SELECT PLOTTED OR TABULATED DATA OPTION. 
2- ENTER THE REQUIRED DATA {SEE ABOVE) WHEN PROS? DEO 
3. TABULAR DATA MAY Be OBTAINED Sean PLOTTING. 


OPTIONS./ GCOARENTS: 


1. TRANSIENT RESPONSE MAY BE OBTAINED FOR IMPULSE, SYEPI eee 
RAMP INPOTS. 

2- A CONTINUOUS QUTPUT IS PLOTTED FOR LAPLACE TPRANSw 
FUNCTIONS. FOR DISCRETE SYSTEMS) fee toe See 
PLOTTED ONLY AT INTERVALS OF THE SAMPLE PERIOD. 

3. Q-4 LINES OF HEADING MAY BE SPECIFIED FOR EACH Pram 
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DISSPLA AND GRAPHICAL OUTPUT DEVICES 


The main purpose of DACSAP is to provide the user with 
graphical systems analysis tools such as root locus plots, 
frequency response flots and time response plots. DACSAP 
takes in user specified parameters and performs calculations 
which result in data which is to be plotted. The actual 
plotting, however, is done by a library of graphics sukrou- 
tines called DISSPLA (Display Integrated Software System and 
Plotting Language). DISSPLA is a web of FORTRAN subroutines 
to which plotting parameters are passed and which, in turn, 
provide the instructions necessary to produce high quality 
plots, graphs, charts or other specialized graphics. 

DACSAP passes data which is to be plotted to one of its 
meeting EtOutaines;  RiLCS, BPLOT, NYPLOT, NICPLT, or TPLOT. 
These routines, then, pass the parameters to DISSPLA by 
Beering a Variety Of subroutines in the proper sequence to 
produce the desired plot. By €Xamining one of the plotting 
routines, the reader can see that, in some cases, no more 
than 30 instructions are required to produce a high quality 
plot with curve smoothing, plotting surface grid, labelled 
axes and headings. 

In general, DISSPLA subroutine names are indicative of 
their function. For example, the subroutine CURVE passes to 
DISSPLA the array names containing the data which defines a 
curve. XNAME and YNAME pass the x and y axis latkels, 
respectively. GRAF, XLOG, and POLAR define the dimensions 
for rectangular, semi-log and polar plotting surfaces, 
respectively. Most of DISSPLA's subroutines follow this 
pattern, thus, making the graphics routines very easy to 


implement. Even if the user is unfamiliar with DISSPLA, he 


i 


could examine one of DACSAP's flotting routines and foliow 
the method used to produce a flot. 

In general, DACSAP controls the size, Shape, scaling and 
nature of information presented on all plots. Thus, the 
user is faced with far fewer questions from the program, but 
because all plotting decisions are made by the program, he 
has little control over the lcok of the plot. The only 
parameter controlled by the user is the plot heading. 
DISSPLA allows up to four lines of heading per plotted page. 
DACSAP accommodates this capability, restricting the length 
of each Jine to 32 characters. 

The orientation of root locus, polar and nagnitude=p pia 
frequency response flots are changed based on whether a 
heading is selected or not. These plots are presented side- 
ways on the TEK 618 when a heading is selected so that the 
long side of the plot is aligned with the tong edge of the 


screen. In this way, hardcopy prints of the screen better 
fill the page. All plots are presented upright if no 
heading is selected. If the user intendS on running a 


multitude of plots, conducting analysis on the screen, it is 
recommended that the plots be done without headings. This 
reduces the plotting time and presents the plots upright. A 
heading may be added at any time and the plot regenerated 
with a heading if a hardcopy is desired. 

DISSPLA produces plots on a number of output devices. 
The user of DACSAP is given four output device options; TEK 
618, TEK 4113/4114, IBM 3278 and DISSPLA metafile. Since 
DACSAP is meant to be an interactive design tool, it is best 
utilized when output is to a high speed, high resolution 
device such as the TEK 618 or TEK 4113/4114, using a high 
speed moden. The TEK 4113/4114 may be too slow for produc- 
tive interactive design if used with baud rates of 1200 or 


less. 
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Output to the IBM 3278 is in the form of a low resolu- 
tion, unlabelled, Hine printer style plot. The general 
shape of curves can be observed, but real analysis of these 
plots is practically impossible. However, this option may 
be useful if the user wishes to input the transfer functions 
for a large systen, make changes to a system or obtain 
jmapular output, and do so without i gain de Graphics 
terminal. 

The option to output to a metafile is available to those 
users who wish to create publication quality plots using a 
postprocessor and high resolution plotting device such as a 
VERSATEC electrostatic plotter. The plots, in this case, 
are sized so as to fit comfortably onan 8.5 by 17 inch 
sheet of paper. Ali figures in this thesis were produced 
using the metafile option. This option may not be used for 
analysis, however, Since outfut does not appear at the 
terminal. Separate terminal sessions are reguired; the 
Tirst to analyze the system and the second to produce plots 


based on parameters obtained in the analysis session. 


oy 


APPENDIX D 
PROGRAMMER*S GUIDE 


This appendix contains detailed information about each 
of DACSAP's subroutines. The main purpose of this appendix 
‘is to provide supplemental information to that found in 
Chapter 3 of this report. Chapter 3 describes the capabili- 
ties and options available from the major routines and 
outlines some of the equations used. This guide outlines 
the specific functions of the driving program, » DACSAP}aygaaa 
each of the subroutines. DACSAP is presented first, then 
the subroutines are presented in alphabetical order. The 
outline for each contains definitions of the main variables 
and arrays, notes about the routine's operation anda list 
of subprograms called by that routine. Table I iS a cross 
reference of subprograms showing the calling routines for 
each subprogran. 

It is hoped that the combination of the program desSeeae 
tion found in Chapter 3 and the more detailed subroutine 
descriptions found in this appendix will help anyone inter- 
ested in interpreting the actual FORTRAN code for the 
program. 


D.1 DACSAP 

_DACSAP is .the_ driving program which presents the main 
option menu and calls subroutines based on the user's selec- 
tion from that menu. 


Definition of variables 


IRS = Problem restart flag . 

ITRM = Graphical output device number 

FD = Printed output device number 1 a ee . 

TINIT = Transfer function farameter initialization flag 
FINIT = Frequency response parameter initialization flag 
IMINIT = Time response parameter initialization flag 
RINIT = Root locus parameter initialization flag 
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TABLE I 


Subroutine Cross Reference 
GALLENG ROUTINES 


dW3GM 
ANIM 
GAAVWSL 
Old: 
aro Td 
lcd i a 
YaAWI AL 
WILEVL 
OUddVaL 
AWILS 
SOTLY 
VLVdihd 
LOIdAN 
iad nN 
d' IH 
yOuds 
NO 

WI LOVd 
LOWOWG 
OUAoVd 
dOOTD 
NESE) 
LOTdd 


dwSowvd 


ATANH 


eeTe Nm | J 


CLINE 


Pai fe ey aie toe | (aa Seam | | 


BPLOT 


sled sa ne ko) | Ag i es) p | =| | 

AOH SES SUM ADAHHA MH MMAH MN AZHMOSHMBAOH Me ab 

ORBDHOVUVUAOHUYOHOVHA KRCZAAHUVUMASMHA MD MOS 4ZaSs 

OmMHbm =MMAZOaNANn AD DHAAHAHHMAHheASaON aH AH 

WH OVO YU a Oe Oe eee ciaeun. ac 

UO dc « Ol fy gH aH PC eee) PaO) CORE te eh 2) El Et 4 E N 
aes xy mea ao mo VY) Et Et Es 4 


SUNT LNOYENS 


a9 


Notes 


1) The routine insures that the transfer function 
parameters have been defined before any analysis 
oe are called by checking the condition of 


2) Analysis routine arameter initialization flags 
(RINIT, FINDT, TMi) are used to determine where, 
within the called subroutine control will be 
transferred. If the flag is 6, control will go to 
the parameter input rcutine. li the flag ~220 
control will go to the parameter change routine. 
Values greater than 1 are temporarily uséd to indi- 
cate special event sequences. 


Subprograms called 


DACFRQ, DACROT, DACTIM, HELP, RNULLI, SCREEN, MTINPUT, 
TSAVE DISSPLA Subroutines. 


D.2 FUNCTION ATANH 


This function evaluates the inverse hyperbolic tangent 
of a real number, -1<x<1. 


Subprograms called 


None 


D.3 SUBROUTINE BFILL 


This subroutine fills a Elock with the closed loop 
transfer function of another systen. 


Definition of variables 

BN = Block number 

NNM = Order of the numeratcr polynomial 
NDN = Order of the denominator polynonial 
BNUM = Array of numerator coefficients » 
BDEN = Array of denominator coefficients 
Notes 


1) See Subroutine CLOOP for a definition of closed 
loop parameters. 


2) This routine is called anytime the _user chooses to 
input a block's transfer function from a data file 
or when expanding the system to an outer loop. 

Subprograms called 


None 
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D-4 SUBROUTINE BPLOT 

tits "SupEOUeinemmecontalns alle@or the calls to, DISSPLA 
pour thes necessary to create a Bode freyuency response 
plot. 


Definition of variables 


= Low end of frequency range 
WMAX = High end of fréquency range 
AMIN = Minimun anee e . 

WFREQ = Array of Frequencies to be Dee Led 

OMAG = Array of magnitudes (dB) to be plotted 
OPHSD = Array of phase values (deg) to be plotted 
LINES = Array of plot heading character strings 


Notes 


1) Subroutine PLTSCL is called prior to this _one to 
calculate plot dimensions based onthe values in 
the arrays QPHSD and QMAG. 


2) The COM ne oes both magnitude and phase on the 
Same plot with_magnitude labelied on the left side 
of the plot and phase labelled on the right. 


Subprograms called 
CLINE, DISSPLA subroutines 


D.5 SUBROUTINE CLIN 


CLIN is the subroutine used to prompt the user for the 
closed loop parameters and is used whenever the program is 
going to save the closed loop transfer function or going to 
calculate the closed loop transfer function of an inner loop 
prior to a system expansion. 


Definition of variables 


KG = Forward path gain 

KH = Feedback path gain a 

TYPFB = Type of feedback (F=positive, N=negative) 

Ii = Flag which determines how user will be prompted 


Notes 

1) This routine checks to see if the closed loop 
parameters. have been reviously defined by some 
Other routine. If so he user 1S simply presented 
with the eee Ty held values and SE the oppor- 
tunity to change them, if so desired. 


2) The subroutine has an integral change and correc- 
tion routine. 


Subprograms calied 
RCHAR, RNULLI, RNULLR, SCREEN 
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D.6 SUBROUTINE CLINE 
This routine a one row of the 2-D array LINES into 
a 1-D array suitable for use by DISSPLA subroutines. 


Definition of variables 


ILN = Line number (row numter of LINES) 
LINES = 2-D aliay sot oe) Ome dene = 
LINE = 1-D array Containing one ron OCfee IES 


Notes 


1) The dollar sign delimiter is placed at the end of 
the character string as reyuired py DISSPLA. 


Subprograms called 


None 


D-7 SUBROUTINE CLOOP 


This subroutine calculates the coefficients of the 
system's closed Loo transfer function. The inputs for the 
routine are the system's open lcop transfer function and the 
individual block transfer functions. 


Definition of variables 


Order of the closed loop numerator 

Order of the closed loop denominator | 

CNM = Array of closed loop numerator coefricients 
CDEN = Array of closed loop denominator coefficients 
CLRIN = eae of closed locp numerator roots 

KN = Closed Loop numerator gain constant 

NNUM = Order of the open loop numerator 

NDEN = Order of the_open loop denominator. 

NM = Array of open loop numerator coefficients 


O) 
tH 
< 
=, 
eal 


DEN = Array of open loop denominator coefficients 
aa = Array of open loop numerator roots 


AGEday 708 ie loop denominator roots 

KG = Forward path gain 

KH = Feedback poe gain ; 

TFID = Array of transfer function I.D. numbers 
{see TINPUT) 
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D. 8 


Notes 
1) The following equation was developed, using Mason's 


ain rite, cor calculating tne closed ioop transfer 
WG LOM oC Lr ¢ 


KG x NG x DH 


CEnr — 
(DG x DH) + (KG x KH x NG x NB) 


KG x NG x DH 
DGEH + (KG x Kad x NGH) 


where NG Forward path numerator 


NH = feedback path numerator 
DG = forward path denominator 
DH = feedback path denominator 
NGH open loop numerator 


DGH = open loop denominator 
2) Multiplication of NG x DH is_ done PY combining 
Enel Eeroots, then calculating the resulting polyno- 
mial coefficients. 
Subprograms called 


MAKPOL 


SUBROUTINE DACFRO 


This subroutine is the frequency response parameter 


input routine. It creates all frompts and menus required to 
interactively input and change frequency response analysis 
parameters. 


Definition of variables 


KG = Forward path gain 

KH = Feedback path gain a 

TYPFB = Type of feedback (F=positive, N=negative) 

K = Open loop gain a: 

TYFQ = ee of frequency (w=real, v=ficticious) 

T = Sample period 

WHIN = Low end of frequency range 

WMAX = High end of frequency range 

WNY = ae frequency (discrete systems) 

NL = Number of lines in the plot heading ; 

LINES = Array of plot heading character Spee ads 

TYPBOD = Type Bode flag (O=open loop, 1=closed loop) 

TYPP = Type of flot eee 

CAC = Change_and correction od same aS parameter 
initialization flag, see DACSAP) 


RN = Return number, determines where control should 
return to after correction routine 

SN = Screen number 

M1 = Change menu flag 


M = Parameter change flag 
MM = Parameter correction flag 
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Notes 


1) Open or closed loop parameters are input based on 
the value of TYPBOD. 


2) Corrections routine 1s used to display to the user 
the parameters which have just been entered. 
Commanded GO TO statements are used to select the 
correct line number ir the input routine when 
corrections are necessary. 


3) Change _routine calls the parameter change menu, 
then selects the corrections. routine if parameters 
are to be changed. Subroutine TINPUT is called if 
changes to transfer functions are desired. 


4) When no further changes _are desired, the calcula- 
tion routine FREOR iS cailed. 


Subprograms called 


FREOR, HELP, CHAR, RLINE, RNULLI, RNULLR, SCREEN, 
TABFRO, TINDOT 


D.9 SUBROUTINE DACRUT 


This subroutine is the root locus parameter input 
routine., It creates all prompts and menus reguired to 
1h te rea cay input and change root locus analysis parame- 

ers. 


Definition of variables 


KMIN = Minimum open loop gain value 
= Maximimum open loop gain value 
TYPFB = Type of feedback fesse N=negative) 
dee C 1c abk fag 
KTICI = Gain increment for tic marks 
T = Sample period (discrete systems) 
XMAX,XMIN = x-axis dimensions 
YMAX,YMIN = y-axis dimensions : 
NL = Number of lines in the plot heading 
LINES = Array of plot eee character strings 
CAC = Change_and correction aoe feats as pabanetes 
initialization flag, see DACSAP) 
SN Screen number 
M1 Change menu flag 
M = Parameter change flag 
MM = Parameter correction flag 
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Hotes 


1) Value _of ITIC determines whether tic marks will be 
placed on the root locus plot. 


2) Corrections routine is used to dispiay to the user 
the parameters which have just been entered. 
Commanded GO TO statements are used to select the 
correct line number inthe input routine when 
corrections are necessary. 


3) Change routine calls the parameter change menu, 
then selects the corrections, routine if parameters 
are to be changed. Subroutine TINPUT is called if 
changes to transfer functions are desired. 


4) When no further changes _are desired, the calcula- 
tion routine RTLCS iS called. 


Subprograms called 
HELP, KCHAR, RLINE, RNULLI, RNULLR, RTLCS, SCREEN, 
TINDOT 


D.10 SUBROUTINE DACTIM 

This subroutine is the time response parameter input 
routine... It creates all prompts and tm§enus reguired to 
eco 1 Vety input and change time response analysis paran- 
eters. 


Definition of variables 


IN = Vee input designator 

UMAX = Max amplitude of infut 

KG = Forward path gain 

KH = Feedback path gain a 

TYPFB = Type of feedbacx (P=positive, N=negative) 

Tt = Sample period (discrete systems) 

TMAX = Length ou time over which response is to be 
analyze 

NL = MunMeee or lines in the plot heading 

PENNE S = AEray Of splot Devan character strings 

CAC = Change and Correction 4 heel aS parameter 
iMimeelalazation Elag, see DACSIA PR) 

RN = Return number, determines where control should 

return to after correction routine 

SN = Screen number 

M1 = Change menu flag 

M = Parameter change flag 

MM = Parameter correction flag 


oO 


Notes 


1) The prompts used depend on whether the system is 
continuous or discrete. 


2) Corrections routine is used to display to the user 
the parameters which have just been entered. 
Commanded GO TO statements are used to select the 
correct line number in the input routine when 
corrections are necessary. 


3) Change routine calls the parameter change menu, 
then selects the corrections routine if parameters 
are to be changed. Subroutine TINPUT is called if 
changes to transfer functions are desired. 


4) When no further changes are desired, the calcula- 
tion routine TIMER is called. 


Subprograms called 


HELP, RCHAR, RLINE, RNULLI, RNUOLLR, SCREEN, TABTIM, 
TIMER, TINPOT 


D.11 SUBROUTINE FCN 

This subroutine defines the_canonical form of the state 
equations, based on the closed loop transfer function, which 
are used to calculate the s-domain time response. 


Definition of variables 


XDOT = Array of values for each differential equation 
at time, ITIME ; | 

ITIME = Time values at which calculations occur 

X = Array of the current values of each state 

MN = Number of states : 

VIN = Input® amplitude at tine, erie 


Notes 


1) This ,routine,is called b the IMSL differential 
eguation solving routine DVERK.- 


Subprograms called 
SINPT 


D.12 SUBROUTINE FOSWCH 

This subroutine performs the bilinear frequency trans- 
Oeil en for frequency response analysis in the w- ofr 
w'-domains. 


Definition of variables 
FOF = Type-of-frequency-transform designator 
NOF = Number of frequency value in the array WFREQ 


WFREQ = Array of fre€cuency values to de plotted 
T = Sample period 
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Notes 


1) The value of FQF .is set by the eee routine, 
FREQR, and determines in which domain and in which 
direction the transform is being performed. 


Subprograms called 


None 


D.13 SUBROUTINE FREOR 


= <> <a ee eee ee oe 


This subroutine is the main calculation routine for 
frequency response analysis. Its purpose is to evaluate the 
Bee apie ate ERawolLem@e SuUnCEION at jo, Wy ob jv’, dagepending on 

e domain. 


Definition of variables 


IYPBOD = Type Bode flag (0=open loop, t=closed loop) 


TYPP Type of rlot flag nee 
NCOF = Array of numerator coerficients 
DCOF = Array of denominator coefficients 


NOF = Number of frequencies in array WFREQ 
WFREQ = Array of frequencies for evaluation 
CNUMR,QNUMI = Real and imaginary value of the 
humerator at jw 
QDENR,QDENI = Real and imaginary value of the 
_ denominator at ae 
QMAG = Magnitude of the transfer function at jw 


QOR,QI = Real and imaginary part of QMAG 
QPHSR = Arrav of phaSe values eee 
QPHSD = of phase values (degrees 


D Array ¢ 
FSF = Phase Shift flag : 
IDOM = Domain of transfer functions 
PM = Phase margin 
GM = Galn Margin. 
EMF Phase wlan gin flag 
GMF Gain margin flag 
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Notes 


1) Subroutine PVAL is _ called _ to evaluate the numerator 
and denominator polynomials at jw, jy or jp". 


2) Complex division is perfomed to determine the 
values of QR and QI. 


3) QR and QI are evaluated to determine phase values. 

4) TYPBOD determines whether NCOF and DCOF are coeffi- 
cients of the open lcop or closed loop transfer 
function. 

5) If the frequency range.chosen by the user allows 
the gain and phase margins are calculated. PMHUF and 
GMF are. set if ain and phase crossovers occur. 


The condition of these flags determines whether PM 
and GM are evaluated. 


6) The value of TYPP determines what type of plot will 
be used to depict the frequency response. 


Subprograms called 
BPLOT, CLOOP, FOSWCH, NICPIT, NYPLOT, PlLisel sees. 


D.14 SUBROUTINE HELP 
This subroutine provides information about the operation 
of each of DACSAP's major subroutines. The routine may be 
accessed from any of the the program's option menus. 
Definition of variables 


I = Calling routine designator 
FD = Output device number 


Notes 


1) The_routine for which information is to be provided 
is determined by the value of I. 


Subprograms called 
SCREEN 


D.15 SUBROUTINE KRANGE 


This subroutine determines the number of roots of a 
po? none which are out of the plotting area of a root 
ocus plot. 

Definition of variables 


N 


Order of the Doe one being evaluated 

Array of coefficients cf the characteristic 
polynomial. ‘ ; 

XMIN, XMAX = x-axis dimensicns of the root locus plot 
YMAX = y-axis max dimension 

Z = Array of complex roots of polynomial, A 

X,Y = Real and imaginary parts of Z 

oft = Number of roots which are outside the 

plotting area 


108 


Notes | 

1) The routine is called Ly  RTLCS to determine if the 
user's ane of gain values is too large for the 
chosen plotting area. 

Subprograms calied 


None 


D.16 SUBROUTINE MAKPOL 


This subroutine calculates the complex coefficients of a 
polynomial given the roots of that polynomial. 


Definition of variables 


N = Order of the eet ae 

R = Array of complex roots of the polynomial 

C = Array or complex coefficients of the polynomial 
Notes 


1) The value of the coefficient of the highest power 
of the unknown is always unity and is not returned 
by the routine. 


Subprograms called 


None 


D.17 SUBROUTINE NICPLT 


This subroutine contains all of the calls to DISSPLA 
subroutines necessary to create a log magnitude-phase plot 
for use in Nichols frequency response analysis. 


Definition of variables 


YORIG,YMX = Dimensions of the phase axis | 
AMIN, AMAX -= Dimensions of the magnitude axis 
opHsDb = Array of freguencies to be Pewene a 

QMAG = Array of magnitude to be plotted ; 
LINES = Array of plot heading character strings 


Notes 

1) Subroutine PLTSCL ais .called prior to this one to 
calculate the lot dimensions based on values in 
the arrays QPHSD and QMAG. 

2) The plot is rotated on its side if a headin 
exists. This will cause the plot to better fil 
the screen and cause any printout of the screen to 
better fill the page. 

Subprograms called 


CLINE 
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D.18 SUBROUTINE NYPLOT 

This subroutine contains all of the calls to DISSPLA 
subroutines necessary to create a polar plot of magnitude 
and phase for use in Nyquist frequency response analySis. 


Definition of variables 


ORMIN,QRMAX = Max and min values of the array QR 
QIMIN,QIMAX = Max and min values of the array QI 

OM = Array of magnitude values 

QR,QI = Real and swag ea Part. Ofer 

RL = Resulting length of axis 

RS TERS Axl Ss St epest 7c ae — 

XDIST @BDIST = Amand 7 pose of Origin of pedares. 
LINES = Array of plot heading character strings 


Hotes 


1) The scaling routine is used to determine the dimen- 
Sions of the plotting area based on the arrays QR 


and Qi. 
2) The plot is rotated on its side if a_ headin: 
exists. This will cause the plot to better fil 


the screen and cause any printout of the screen to 
better fill the page. 


Subprograms called 
CLINE 


D.19 SUBROUTINE PLTSCL 


This subroutine calculates the dimensions for pe 
response plots based on the arrays of values to be plotted. 


Definition of variables 


AMIN, AMAX = Min and max values of the array QMAG 
PMIN, PMAX = Min and max vaiues of the array QPHSD 
YORIG, YMX = Min and max values of the phase axis 


QMAG = Array of magnitude values (dB) 
QPHSD = Array of phase values (degrees) 


Notes 

1) The routine insures that numbers that are_ easy to 
work with are used on the axes of frequency 
response plots. 

Subprograms called 


None 


D. 20 


_This subroutine calculates the open loop gain, damping 
Bat2O and natural Erequency of any point on a root locus 


DLot. 


Definition of variables 


kone POU ooo by the user 
T = Sample perio 

SIG = Damping 

WD = Damped natural frequency 


aes = Domain of the transfer functions 


Damping ratio 

W Undamped natural frequency 

K Open itoop gain 

XD, YD = Comparison x and y values 

TOL = Tolerance value for reading X and Y from the 
root locus plot 


Notes 


1) The routine solves explicitly for values in the s-, 
w- or w'-domains. 


2) w- andw'-plane natural freyuency, YW, can be 
expressed as 


Wee Sot wey te 


3) w- and w'-plane damping ratio, D, can be expressed 


as 
D= { SIG /W | 
4) w- and w'-plane damping and damped natural 
frequency are 
w-plane w*'-plane 
SIG = aA tanh-i X SIG = An tania ( 6 172) 
WD = (Z/T) tan-!? Y WD = (Z/T) tan-!t (Y T/2) 


5) The z-plane iteration routine calculates X and Y 
values for all vaiues cf damping ratio and natural 
frequency compares them to the user's X and Y 
value an aaSP ays the values of W and D which 
correspond to the closest match. 


Subprograms called 
None 


D.21 SUBROUTINE PVAL 


This subroutine evaluates a polynomial at any complex 
value of the variable. 


Definition of variables 


A = Array of polynomial coefficients 

NN = Order of the polynomial 

S = Complex value at which the polynomial is being 
evaluated 

P = Complex value of the pclynomial at S 

PR,PI Real and imaginary parts of S§S 

VR, Vi Real and imaginary parts of P 


Subprograms called 


None 


D.22 SUBROUTINE RCHAR 


; This subroutine 1s used to read a literal reply to an 
interactive prompt allowing for a recovery from an inadver- 
tent null string entry. 

Definition of variables 


IANS = Character response to the prompt 
COUNT = EXKror count 


Notes 


1) Two null string entries in a row and program execu- 
tion will end. 


Subprograms cailed 


None 


D.23 SUBROUTINE RLINE 

This subroutine reads in a character string from the 
console allowing fora recovery from an inadveftent null 
StLiNG@ eCnvury. 

Definition of variables 


ILN = Number of lines in the array LINES 
LINES = 2-D array Of Chakdcters 


Notes 


1) Two null string entries in a row and program execu- 
tion will end. : 


subprograms cailed 


None 


D.24 SUBROUTINE RNULID 


Titsesupsoutineeresused tO (readva double precision real 
number reply to. an interactive prompt ailowing for a 
recovery from an inadvertent nuli string entry. 

Definition of variables 


ANSR = Double precision response to the prompt 
COUNT = Error count 


Notes 


1) Two null string entries in a row and program execu- 
tion will end. 


Subprograms called 


None 


D.25 SUBROUTINE RNULII 


; This subroutine is used to read an integer reply to an 
interactive prompt allowing for a recovery from an inadver- 
tent nuil string entry. 

Definition of variables 


IANS = Integer pes ons to the prompt 
COUNT = Error coun 


Notes 


1) Two null string entries in a row and program execu- 
tion will end. 


Subprograms called 


None 


D.26 SUBROUTINE RNULIR 


This subroutine is used to read a Singie precision real 
number reply to, an interactive prompt’ allowing for a 
recovery from an inadvertent null string entry. 


Definition of variables 


ANSR = Single precision response to the prompt 
COUNT = Error count 
Notes 


1) Two null string entries in a row and program execu- 
tion will end. 


Subprograms called 


None 


De 


This subroutine is the calculation and plotting routine 


for root locus analysis. 


Definition of variables 


FOUT = Number of zeroes outside the ep ens area 
KMX = Calculated max value of gain which will be used 
tor plotting purposes 


KLIM = Interim value of KMX 

KOUNT = Number of roots plectted 

NNUM = Order of the open lcop numerator 

NDEN = Order of the open loop denominator. 

NM = Array of open loop numerator coefficients 


NUM = Array of open loop numerator coefficients with 
leading zeroes — 

DEN = Array of open nook dé€nominator coefficients 

Z = Array of complex por ya cw roots (output of ZPOER) 

XPOLE, YPOLE = Real and imaginary parts o£ system poles 

XZERO,YZERO = Real and imaginary parts of system zeroes 

LINES = Array of plot heading character strings 

XMIN, XMAX = x-axis plot dimensions 

YMIN,YMAX = y-axis plot dimensions . ar 

A = Array of characteristic polynomial coefficients 

K = Open loop gain_ used in calculating A_. 

KP = Particular value of gain to be highlighted on plot 

KTMIN,KTMAX = Range of gain values for tabulated output 

KTINC = Gain increment for tabulated ouput 


Notes 

1) This routine is a calculation, plotting and tabular 
data routine. This is done only for _ root_locus 
analysis because roots are calculated for only one 
set of roots at a time, lotted or tabulated and, 


then, gain is incremented and the process begins 
again. Arrays containing all of the roots and gain 
values are hot formed because they would take up 
too much memory. 


2) The plotting routine contains all of the DISSPirA 
subroutine “calls necessary to draw. a root iocus 


lot except for the discrete domain plots. In 
hese caseS, TEMPLT cr WTEMP are also called to 
draw the ines of constant damping ratio and 


natural frequency. 


3) The beginning of the routine adjusts the max gain 
value 1f the user has chosen one which is so large 
that the roots cannot be plotted in the plotting 
area selected. 


Subprograms called 


Fae ae KRANGE, PIDATA, 
a 


RCHAR, RNULLI, RNULLR, SCREEN, 
TEMPLT, WTEMP, ZEOLR (IM 


CH 
SL) 
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D.28 SUBROUTINE SCREEN 


. This subroutine contains the instructions required to 
display, at. the terminal, alli option menus and parameter 
change and display screens. 


Definition of variables 

SN = Screen number 

Notes 

1) A description of. variables used in, each of_ the 
screens may be found in the associated calling 
routine. 

Subprograms called 


None 


D.29 SUBROUTINE SINPT 
This subroutine calculates the magnitude of the input at 
any time for time response calculations. 


Definition of variables 


= Type-of-input designator | 
ATTIRE = Time at which calculation is taking place 
Time interval for ITIME . 
Max time (sec) over which response is observed 
Max amplitude of the input. 
UIN = Amplitude of the input at time ITIME 


Notes 


1) Variable IN is established by DACTIM and determines 
the type of input (stef, impulse or ramp). 


+ 
= 
od 
< 
Welw 


Subprograms called 


None 


D.30 SUBROUTINE STIME 

This subroutine calculates the coefficients of the 
canonical form of the state, equations based on the closed 
loop transfer functicn coefficiénts. 


Definition of variables 


TOL = Integration tolerance 

TIME = Array of time values : 

ITIME = Time at which calculation is taking place 
X = Array of state variables 

CLNN = Order of the closed loop numerator 

CLND = Order of the closed loop denominator 


CNM = Array of closed loop numerator coefficients 
CDEN = Array of closed loop denominator coefficients 
CNUM = Array of closed loop numerator coefficients 


with leading zeroes 


A = Coefficients of the companion matrix 
B = Coefficients of the output equation 
AMP = Array of output amplitude values 


1S 


Notes 


1) The array A contains the coefficients of the bottom 
row of the canonical form of the state matrix: 
These coefficients are based on the denominator 
coefficients of the clcsed loop transfer function. 


2) The remainder of the coefficients in the comfanion 
matrix are i's in the superdiagonal and are taken 
care of in the subroutine FCN. 


3) The array B contains the coefficients of the output 
equation and is based cn the numerator coefficients 
of the closed loop transfer function. 


4) The last part of the routine contains the _loop 
which calls the IMSL differential equation solving 
routine DVERK. | DVERK solves the set of sinulta- 
neous differential eguations described in FCN, 
which uses the coefficients calculated in STIME. 


Subprograms called 
DVERK (IMSL), FCN (through DVERK) 


D.31 SUBROUTINE TABFRQ 

This subroutine is the frequency response tabulated data 
parameter input routine. It _ creates all of the prompts and 
menus required to interactively input and change the parane- 
ters used to tabulate frequency response data. 


Definition of variables 


K = ORE loop gain. 
KG orward path gain 

i Feedpacw patn cain a. 
TYP = Type of feedback (N=negative, P=positive) 
WMIN,WMAX = Range ot Erequencies for plotting 
FMIN, FHAX = Range of frequencies to be tabduiated 
NOF = Number of frequencies to_be tabulated 
NLT = Number of lines for tabulation heading 
THEAD = Array of heading character strings 
IM = Correction routine flag 
FD = Output device number 


Notes 


1) The variables TK, TKG Paton a str yPrb es are local 
variables used to hold. the associated Proaung 
parameters constant durin tabulation. Dees 
ae Keen Gee KHoand S2yYer. fe eng common vari- 
ables) noe be used to sen DMeoetatmon tO ~tne 
routine SCREEN. 


Zz). ere ee ee warameters have previously 
been defined, hen reentry 1s not. required. 
Current parameters are displayed and may be 
changed, if desired. 


3) The correcticns routine is used to display to the 
user, the parameters just entered. Commanded GO TO 
Steele nesmedkhe. USCOM tO) SemeCEm tne COLrect line 
humber in the input routine when corrections are 
necessary. 


Subprograms called 
PRO Chee Uc ei oernnULG AG oCh EEN REAR 


D.32 SUBROUTINE TABTIM 

This subroutine is the time response tabulated data 
Boeametecr input routine.  ~@t creates ail of the prompts and 
nenus eeu ed to interactively input and change the parame- 
ters used to tabulate time respense data. 


Definition of variables 


IN = Ce een Ste designator 

UMAX = Max amplitude of the input 

KG = Forward path gain 

KH = Feedback path gain Pe. 

TYPFB = Type ort feedback (N=negative, P=positive) 

T = pou eriod (discrete systems) 

TMAX = Aax time over which the response is observed 

NLT = Number of lines for tabuiation heading 

THEAD = Array of heading character strings 

UIN = Amplitude of the input at_time TIM 

TINT = Interval between tine value in TIME 

TIM = Current value of time in TIME 

TINE = Array of time values used in calculations 

AMP = Array of output amplitude values for each 
Erneerne LLM . 

at @Conpeecctton routine flag 

=D Output device number 


Notes 


1) 


2) 


3) 


The variables TTIN, TTIUMAX, MTITKG ren TT Ye Re 

TTT,  “anGel TIMNAX Jamemlocas variables used to hold 
the associated piotting parameters constant during 
tabulation. In this way, IN, UHAX, KG, KH, 2 voi 

T and TMAX (existing ccmmon variables) may be used 
to send information to the routine SCREEN. 


If time response parameters have ee been 
defined, then reentry is not reguired. urrent 
aneneTers are displayed and may be changed, if 
esire e i 


The corrections routine is used to display to the 
user, the parameters just entered. ~Comanded GO TO 
Statements are used to select tne correct line 
Nhumber in the input routine when corrections are 
necessary. 


Subprograms called 
FREOR, RNULLI, RNULLR, SCREEN, Shean 


D.33 SUBROUTINE TEMPIT 


Thies 


subroutine calculates and plots lines of constant 


damping ratio and natural frequency for z-plane root locus 


ELGES. 


Definition of variables 


XAF = Length of root locus x-axis ; 

YAF = 10% of length of root locus y-axis _. 

ISTEP = Step factor for interval between lines of 
constant natura) 2 tee 

JSTEP = Step factor for interval between lines of 
constant damping ratio_- 

P = Array of label values for lines of constant 


WW = 


natural frequency 
Natural frequency 
First line of natural frequency to be labelled 


s-plane coordinate corresponding to some value 


of natural frequency and damping ratio 


RT,IT = Real and imaginary parts of 5 


x, ¥ 
D= 


Z-EbontosltOrMeon as 

= PEL of real and imaginary parts of Z 
(plo LR rer 

Dap Ea tae 
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Notes 


1) 


2) 


3) 


The value of S associated "with givei values ofr 
natural frequency and damping ratio is 


5S. = suet ($2 aoe ec 


The associated value of Z 1S 
Z= exp (sf) 


where T is the sample period. The complex values 
of Z are calculated and plotted by this routine. 


The density of lines lotted varies depending on 
the dimensions of the otting area. The instruc- 
tions which establish values for ISTEP and JSTEP 
determine this density. 


The array P is used so that each line is labelleda 
max of one time. When a line is labelled, the 
aoe ape tag clement in P is set to zero which, 
thereaiter, indicates that that line should not bs 
labelled again. 


Subprograms called 


None 


D-34 SUBROUTINE TIMER 


TIMER 
fron Of t 


Defi 


AMP 
AMPH 
1g 


is a. housekeeping routine called prior to execu- 
he main time response calculation routines. 


nition of variables 


Domain of transfer functions 

Time interva J : 

Max time over which time response 1S observed 

Number of points to be plotted 

Number of periods between O and TMAZ 
discrete Seems) 

Indicates number of periods will be plotted 
{discrete systems) 

= Array of amplitude to be plotted 

IN, ANPMAX = Min and max values in array AMP 
Plot or tabulation OU SO iia Lag 


Notes 


1) 


2) 


The function of this reutine is ads follows: 
a) call CLOOP to caiculate loop transier Luncewmee 


b) call WINV to wand w' transfer functions to z 
1f necessary 


c) establish the interval of time response calcu- 
lations based on the domain 


d) call STIME for continuous systems or ZTIME for 
discrete systems 


e) determine dimensions of plot amplitude axis 

£f) cali the plotting routine TPLOT 

A max of 500 points will be plotted for. discrete 
systems. if more than 500 sample periods are 


covered by TMAK, N1 is increased from 1 as neces- 
Sary and only every Ni sample period is plotted. 


Subprograms called 
CLOOP, STINE, TPLOT, WINV, ZTIME 
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D.35 SUBROUTINE TINPUT 

This is the transfer function input and manipulation 
Lowe ane. It creates all of the BEomp S and menus required 
to interactively input and chanye block diagram transfer 
Dierlton PDabameotersom, the Foutine also contains the instruc- 
tions which multiply together the block transfer functions 
to form the systems! open loop transfer function. 


Definition of variables 


NNM = SE COREG the order of the numerator of 
eac Gis 
NDN = Array containing the order of the denominator 


of each block - 

CONN = Array containing the gain constant of the 
numerators of each Llock . 

COND = Array.containing the gain constant of the 
denominator of each block ‘ 

BNUM = 2-D array of numerator coeffecients for each 


block 
BDEN = oe aura! of denominator coefficients for each 
oc 
RTN = 2-D array of numeratcr roots for each block 
RTD = 2-D array of denominator roots for each block 
RR, R1 = Real and imaginary part of a root (input 


variable) 
RN = Array of the numerator roots of the open loop 
transfer function 


RD = Array of the denominator roots of the open loop 
transfer function —_ 

NM = Array of numerator coefficients of tie open loop 
ehaoecrmrunct 120n = 

NUM = Array of numerator ccefficients of the open 


loop transfer functicn with leading zeroes 
DEN = Array of denominator coefficients of the open 
Toopm@transfer functicn 
NNUM = Order of the open Jcop numerator 
NDEN = Order of the open loop denominator 
A = Ne ey of polynomial coefficients evaluated by 
Z = Array of complex roots returned PO aera 
NORD = ee ae or denominator operation flag (see 
notes 
LPNO = Loop number 
BN = Block number 
NBLKS = Number of blocks J 
IDOM =_Domain of transfer functions _ 
IBF = Intermediate block operation flag (see notes) 
TPORN = Form o£ transfer funetion 
C = Array of coefficients returned by MAKPOL 
CC = pee ey C in reverse order of elements 
DIFF = Difference between NNUM and NDEN 
CAC = eae on corrections flag (same as TINIT, see 
EBN = Expansion block number (see notes) 
CcC = Type of polynomial flag (see notes) 
M = Parameter change flag 
MM = Parameter correction ed crn fe. 
IMINIT, RINIT, FINIT = Parameter initialization flags 
_ (see DACSAP) 
TFID = Transfer function I.D. number 
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Notes 


1) 


2) 


3) 


4) 


2) 


) 


cp) 


8) 


*~a) 


The block transfer fume. ton a routine, Stazuima 
at line 1010, is broken into three parts: 


oe to line 1025 defines the characteristics of 
the block» i.e, the path in which it exists, 
the order of the numerator and de nominator 
polynomials and the form of the polynomials. 


b) From 1025 to 1017. is the input routine cae 
transfer functions in coefficient form. This 
part of the routine also calculates the roots 
of the blocks input in coefficient forn. 


c) Line 1017 to 1154 input transfer functions in 
factored forn. 


Line numbers 16XX deal with changing transfer func- 
t20NS 12n another Loop. The _current loop is saved 
ae the other loop 1s loaded with these instruc- 
10ns. 7 


Line numbers 17XX are used to expand to an outer 
loop. These instructions save the_current loop, 
calls the routine to fcrm the closed loop transfer 
functions and places that transfer function in the 
roper place _in.the outer loop. Then the outer 
loop is loaded if it is in a.data file or Conteen 
is sent to the input rcutine if it is to be entered 
from the console. 


IBF is set by the loop. change routines and used Dy 
the input routine to determine which instructronms 
of the input routine should be executed. 


EBN describes the block in the outer loop where tne 
inner loop closed loor transfer function is_ to be 
placed. This parameter is stored when a loop is 
Saved. When an outer loop is loaded, the closed 
loop transfer function of the inner loop_ is auto- 
matically loaded into the current block designated 


by EBN. 


NORD determines which input instructions should be 

executed when changes are pera made, i.e. numer- 

aoe instructions, denominator instructions or 
Gun. . 


ccC is sent to the SCREEN routine to establish the 
type of polynomial parameters which should be 
displayed. 


ccc = 1 humerator in coefficient form 
CCC = 2 denominator in coefficient forn 
CCCs humerator in factored form 
CCC = 4 denominator in factored form 
The array TFID contains a number for each block 
which identifies the block. 
TFID = 36 coefficient form in forward path 
TFID = 1 coefficient form in feedback path 
TFID = 10 factored form in forward path 
TFID = 11 factored form in feedback path 


Subprograms called 
BFILL, CLIN, CI100P, HELE, MAKPOL, RCHAR, RNULLD, 


VZ2 


RNULLI, RNUULR, SCREEN, ILOAD, TSAVE 


D.36 SUBR 


This subroutine reads data from a data file which has 
been formed by the subroutine TSAVE. 


IO 
| 
3 
bed 

(1% 


Definition of variables 


II = eee ee designator 
LPN = Loop number 


Notes 


1) II determines how much of the data file is read in. 


It = 0 Block and system sections are loaded 
II = 1 Block section es 1s loaded 
II = 2 System section only 1s loaded. Biock 


section is loaded into a dummy array. 


2) The remainder of the variables used in this routine 
are defined in TINPUT {section D. 35). 


Subprograms called 
RCHAR 


Do.37 SUBROUTINE TPLOT 


fils SUDbeuUtane “contains all of the calls to DISSPLA 
subroutines necessary to create a time response plot. 


Definition of variables - 


TMAX = Time over which _ the response is to be observed 
AMPMIN,AMPMAX = Min and max values of the array AMP 
TIME = Array of time values to be plotte 

AMP = Array of amplitude values to be piotted 

LINES = Array of plot heading character strings 


Notes 
1) The plot _is rotated on its side if a headin 
exists. This will cause the plot to better fil 


the screen and cause any printout of the screen to 
better fill the page. 


Subprograms called 
CLIN 


D.38 SUBROUTINE TSAVE 


This subroutine stores system parameters to a data file 
under a user specified filename. | The routine is called by 
the user via he DACSAP main oftion menu or by the program 
me CE NP UT. 

Definition of variables 


II = Block-section-load designator 
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Notes 


1) Data is saved into two sections. The rirst section 
contains either the ofen or closed loop transfer 
function of the system. The second contai noma 
other information concerning the system, inciuding 
the transfer functions of each Dlock. 


2) Data in the first section is recalled when indi- 
vidual blocks are loaded from a data file. The 
Meee 1s recalled when the entire systen 
is loaded. 


3) The remainder of the variables used in this routine 
are defined in TINPUT (Section D.35). 


Subprograms called 
CLIN, CLOOP, RCHAR 


D.39 SUBROUTINE WINV 


This subroutine performs the inverse transform of w and 
w' transfer functions to a z-domain transfer function. 


Definition of variables 


CN = Numerator gain constant 
CL = Denominator gain constant eet 
ae Interim array of numerator coefficients 


Interin aes of denominator coefficients 
CLNN = Order of the closed loop numerator 

CLND = Order of the closed loop denominator . 

CNM = Array of closed loop numerator coefficients 
CDEN = eS of closed Os denominator coefficients 
Z Array of roots from sukroutine ZPOLR (IMSL) 


R = Array of transformed numerator roots 
T = Sample period 
Notes 


1) Given the w (or w') transfer function, G(w), 
(w+ a,) (Ww + ay) .-. 
(w + b,) (w + bi) (w+ b,) .-- 


the following relation was developed for the inverse 


Cw 


transforn 
m-n 
ae yes 1) Il ( (a, $262 + eae 

HC); + C)z + (bj + C)} 
j= 

where C= 7 for w to z transforms 

C = 2/T for wt to z transforms 
2) This routine transforms only the closed loop 


transfer function in Pace rates for time response 
calculations (see ZTIME). ; 


Subprograms called 
MAKPOL, ZPOLR (IMSL) 
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D.40- SUBROUTINE WTEMP 

This subroutine calculates and plots lines of constant 
damping ratio and natural frequency ror w- and wt-plane root 
Wacus plots. 


Definition of variables 


IDOM = Domain of the transfer functions (w or w') 

T = Sample period 

XINT = X-interval between lines of constant natural 
frequency 

XMIN = Min value of root lceccus x-axis 

X = Array of x plotting points 

Y = oa of y eee points |. 

XPINT = Interval between values in the array X 

WN = Natural frequenc . 

YM = YMAX or ABS(YMIN) whichever is larger 

Z = Array of damping ratio values to be plotted 


Notes 
1) For a given x value and natural freguency 
y 2/T tan ([T/2 {WN2 - [2/T tanh-! (xT/2) ]2} V2] 
2) For a given y value and damping ratio 
Pau L/ 2) ; 


Xx = 2/T tanh 
(1 - Z2) v2 


3) Five lines of constant Natural frequency are 
plotted at even intervals over the negative x-axis. 


4) Lines of constant damping ratio are plotted for Z = 
tee - & 2 ee 4 = 


Subprograms called 
ATANH 


D.41 SUBROUTINE ZTIME 


This subroutine calculates the time response of a 
z-domain transfer function to a step, impulse or ramp input. 


Definition of variables 


CLNN = Order of the closed loop numerator 

CLND = Order of the closed loop denominator | 

CNM = Array of closed loop numerator coefficients 
CDEN = Array of closed loop denominator coefficients 
ZNNM = CLNN plus the order of the input numerator 
ZNDN = CLND plus the order of the input denominator 
QNN = Array of numerator ccefficients of the CLTF 


ius input ae 
QDN = Array of denominator coefficients of the CLTF 
Pees dn put on. 
A = Local variable used in the long division calculation 
AMP = Array of amplitude values to be plotted 
TIME = Array of time values to be plotted 


Za 


Notes 
1) The closed,loop transfer function lS combined mya 
ae Zea s Foam of either an impulse, step or ramp 
unction. 


2) Long division is used to calculate the output at 
lntervals of the sample period. 


Subprograms called 


None 
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